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Abstract 

We review the main properties, demographics and applications of binary and millisecond 
radio pulsars. Our knowledge of these exciting objects has greatly increased in recent years, 
mainly due to successful surveys which have brought the known pulsar population to over 
1700. There are now 80 binary and millisecond pulsars associated with the disk of our Galaxy, 
and a further 103 pulsars in 24 of the Galactic globular clusters. Recent highlights have been 
the discovery of the first ever double pulsar system and a recent flurry of discoveries in globular 
clusters, in particular Terzan 5. 
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1 Introduction and Overview 

Pulsars - rapidly rotating highly magnetised neutron stars - have resulted in many applications in 
physics and astronomy. Striking examples include the confirmation of the existence of gravitational 
radiation |314| as predicted by general relativity |312l 131 J] . the first detection of an extra-solar 
planetary system |346l 1244] and the discovery of the first double-pulsar binary system |44l I198| . 
The diverse zoo of radio pulsars currently known is summarized in Figure Q 



Figure 1: Venn diagram showing the numbers and locations of the various types of radio pulsars 
known as of January 2005. The large and small Magellanic clouds are denoted by LMC and SMC. 

Pulsar research has proceeded at a rapid pace since the first two versions of this article |179II18()] . 
Surveys mostly with the Parkes radio telescope |259| . but also at Green Bank |232| . Arecibo 228 
and the Giant Metre Wave Radio Telescope |231| have more than doubled the number of pulsars 
known back in 1997. The most exciting new results and discoveries from these searches are discussed 
in this updated review. 

We begin in Section[21with an overview of the pulsar phenomenon, the key observed population 
properties, the origin and evolution of pulsars and the main search strategies. In Section [3J we 
review present understanding in pulsar demography, discussing selection effects and their correction 
techniques. This leads to empirical estimates of the total number of normal and millisecond pulsars 
(see Section I3~3l and relativistic binaries (see Section 1^4^) in the Galaxy and has implications for 
the detection of gravitational radiation from coalescing neutron star binaries which these systems 
are the progenitors of. Our review of pulsar timing in Section 01 covers the basic techniques (see 
Section l4~2")l . timing stability (see Section rO|l . binary pulsars (see Section l4~4f . and using pulsars 
as sensitive detectors of long-period gravitational waves (see Section r4.5f) . We conclude with a 
brief outlook to the future in Section [5] Up-to-date tables of parameters of binary and millisecond 
pulsars are included in Appendix IA1 
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2 Pulsar Phenomenology 



Many of the basic observational facts about radio pulsars were established shortly after their 
discovery |119| in 1967. Although there are still many open questions, the basic model has long 
been established beyond all reasonable doubt, i.e. pulsars are rapidly rotating, highly magnetised 
neutron stars formed during the supernova explosions of massive (> 5-10 Mq) stars. In the 
following, we discuss the observational properties that are most relevant to this review. 

2.1 The lighthouse model 

Figure [3 shows an animation depicting the rotating neutron star or "lighthouse" model. As the 
neutron star spins, charged particles are accelerated out along magnetic field lines in the mag- 
netosphere (depicted by the light blue cones). The accelerating particles emit electromagnetic 
radiation, most readily detected at radio frequencies as a sequence of observed pulses produced as 
the magnetic axis (and hence the radiation beam) crosses the observer's line of sight each rotation. 
The repetition period of the pulses is therefore simply the rotation period of the neutron star. 
The moving "tracker ball" on the pulse profile in the animation shows the relationship between 
observed intensity and rotational phase of the neutron star. 




Figure 2: GIF movie showing the rotating neutron star (or "lighthouse") model for pulsar emission. 
Animation designed by Michael Kramer. 

Neutron stars are essentially large celestial flywheels with moments of inertia ~ 10 38 kg m 2 . 
The rotating neutron star model [2421 111)5) predicts a gradual slowdown and hence an increase 
in the pulse period as the outgoing radiation carries away rotational kinetic energy. This model 
became universally accepted when a period increase of 36.5 ns per day was measured for the pulsar 
in the Crab nebula |27(J| . which implied that a rotating neutron star with a large magnetic field 
must be the dominant energy supply for the nebula 106 . 

2.2 Pulse periods and slowdown rates 

The present public-domain catalogue |208) . available on-line contains up-to-date parameters 
for over 1500 pulsars. Most of these are "normal" with pulse periods P ~ 0.5 s which increase 
secularly at rates P ~ 10~ 15 s/s. A growing fraction are "millisecond pulsars", with 1.5 ms < 
P << 30 ms and P < 10 -19 s/s. For many years, the most rapidly rotating neutron star known 
was the original millisecond pulsar B1937+21 with P = 1.5578 ms. However, one of the 
recent discoveries in the globular cluster Terzan 5 appears to be rotating even more rapidly, with 
P = 1.4 ms |233) . Confirmation of this exciting result should be published in due course. While 
the hunt for "sub-millisecond pulsars" continues, and most neutron star equations of state allow 
shorter periods, it has been suggested j^ElED] that the lack of pulsars with P < 1.5 ms is caused 
by gravitational wave emission from R-mode instabilities 0. 

As can be seen from the ll P-P diagram" in Figure^ normal and millisecond pulsars are distinct 
populations. The differences in P and P imply fundamentally different magnetic field strengths 
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and ages. Treating the pulsar as a rotating magnetic dipole, one may show |185) that the surface 
magnetic field strength B cx (PP) 1 / 2 and the characteristic age r c = P/(2P). 

Lines of constant B and r c are drawn on Figure^ from which we infer typical values of 10 12 G 
and 10 7 yr for the normal pulsars and 10 8 G and 10 9 yr for the millisecond pulsars. For the rate 
of loss of kinetic energy, sometimes called the spin-down luminosity, we have E cx P/P 3 . The lines 
of constant E shown on Figure show that the most energetic objects are the very young normal 
pulsars and the most rapidly spinning millisecond pulsars. 




10 3 0.01 0.1 1 10 

Period (s) 

Figure 3: The P-P diagram showing the current sample of radio pulsars. Binary pulsars are 
highlighted by open circles. Theoretical models \&4Jj do not predict radio emission outside the dark 
blue region. Figure provided by Michael Kramer. 



2.3 Pulse profiles 

Pulsars are weak radio sources. Measured intensities, usually quoted in the literature for a radio 
frequency of 400 MHz, vary between 0.1 mJy and 5 Jy (1 Jy = 10~ 26 W m~ 2 Hz -1 ). As a result, 
even with a large radio telescope, the coherent addition of many hundreds or even thousands 
of pulses is usually required in order to produce a detectable "integrated profile". Remarkably, 
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although the individual pulses vary dramatically, the integrated profile at any particular observing 
frequency is very stable and can be thought of as a fingerprint of the neutron star's emission beam. 
Profile stability is of key importance in pulsar timing measurements discussed in Section 0] 

The selection of integrated profiles in Figure 0] shows a rich diversity in morphology including 
two examples of "interpulses" - a secondary pulse separated by about 180 degrees from the main 
pulse. The most natural interpretation for this phenomenon is that the two pulses originate from 
opposite magnetic poles of the neutron star (see however |2()9| ). Since this is an unlikely viewing 
angle we would expect interpulses to be a rare phenomenon. Indeed, the fraction of known pulsars 
in which interpulses are observed in their pulse profiles is only a few percent |159| . 



(a) B0031-07 (b) B0329+54 (c) J0437-4715 




(d) B0525+21 (e) J1012+5307 (f) B1831-04 




Figure 4: A variety of integrated pulse profiles taken from the available literature. References: 
Panels a, b, d, f \Wi\l , Panel c ]23j/. Panels e, g, i Panel h ]29lj . Each profile represents 360 

degrees of rotational phase. These profiles are freely available from an on-line database [327]. 

Two contrasting phcnomcnological models to explain the observed pulse shapes are shown in 
Figure El The "core and cone" model |261) depicts the beam as a core surrounded by a series of 
nested cones. Alternatively, the "patchy beam" model 202 109 has the beam populated by a 
series of randomly-distributed emitting regions. Further work in this area is necessary to improve 
our understanding of the shape and evolution of pulsar beams and fraction of sky they cover. This 
is of key importance to the results of population studies reviewed in Section 13.21 

2.4 Interstellar dispersion and the pulsar distance scale 

From the sky distribution shown in Figure El it is immediately apparent that pulsars are strongly 
concentrated along the Galactic plane. This indicates that pulsars populate the disk of our Galaxy. 
Quantitative estimates of the distance to each pulsar can be made from the measurement of 
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pulse dispersion - the delay in pulse arrival times across a finite bandwidth. Dispersion occurs 
because the group velocity of the pulsed radiation through the ionised component of the interstellar 
medium is frequency dependent: Pulses emitted at lower radio frequencies travel slower through 
the interstellar medium, arriving later than those emitted at higher frequencies. The delay At in 
arrival times between a high frequency and a low frequency v\ Q pulse is given |185| by 



At = 4.15 ms x 



where the dispersion measure is 



GHz 
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is the integrated free electron column density n c out to the pulsar at a distance d. This equation 
may be solved for d given a measurement of DM and a model of the free electron distribution. 

In practice, the electron density model is calibrated from the 100 or so pulsars with independent 
distance estimates and measurements of scattering for lines of sight towards various Galactic and 
extragalactic sources [3101 IHH5J . The most recent model of this kind [23 E2| provides distance 
estimates with an average uncertainty of ~ 30%. 



2.5 Pulsars in binary systems 

As can be inferred from Figure ^ only 4% of all known pulsars in the Galactic disk are members of 
binary systems. Timing measurements (see Section^J place useful constraints on the masses of the 
companions which, supplemented by observations at other wavelengths, tell us a great deal about 
their nature. The present sample of orbiting companions are either white dwarfs, main sequence 
stars or other neutron stars. Two notable hybrid systems are the "planet pulsars" B1257+12 
and B1620— 26. PSR B1257+12 is a 6.2-ms pulsar accompanied by at least three terrestrial-mass 
bodies [3461 12441 1345) while B1620— 26, an 11-ms pulsar in the globular cluster M4, is part of a 
triple system with a 1-2 Mj up it e r planet [3171 1171 13151 1284| orbiting a neutron star-white dwarf 
stellar binary system. The current limits from pulsar timing favour a roughly 45-yr low-eccentricity 
(e ~ 0.16) orbit with semi-major axis ~ 25 AU. Despite several tentative claims over the years, 
no other convincing cases for planetary companions to pulsars exist. Orbiting companions are 
much more common around millisecond pulsars (~ 80% of the observed sample) than around the 
normal pulsars (< 1%). Binary systems below the line have low-mass companions (< O.7M - 
predominantly white dwarfs) and essentially circular orbits: 10" 5 < e < 0.01. Binary pulsars with 
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Figure 6: The sky distribution of pulsars in Galactic coordinates. The plane of the Galaxy is the 
central horizontal line. The Galactic centre is the midpoint of this line. Millisecond pulsars are 
indicated in red. Binary pulsars are highlighted by the open circles. The more isotropic distribution 
of the millisecond and binary pulsars reflects the differences in detecting them out to large distances 
cf. the normal population (see Section\^). 

high-mass companions (> 1 Mq - massive white dwarfs, other neutron stars or main sequence 
stars) tend to have more eccentric orbits, 0.15 > e > 0.9. 

2.6 Evolution of normal and millisecond pulsars 

A simplified version of the presently favoured model |23 1951 12851 to explain the formation of 
the various types of systems observed is shown in Figure [7] Starting with a binary star system, 
a neutron star is formed during the supernova explosion of the initially more massive star. From 
the virial theorem, in the absence of any other factors, the binary will be disrupted if more than 
half the total pre-supernova mass is ejected from the system during the (assumed symmetric) 
explosion 120 33 . In practice, the fraction of surviving binaries is also affected by the magnitude 
and direction of any impulsive "kick" velocity the neutron star receives at birth from a slightly 
asymmetric explosion |12()l 120) . Binaries that disrupt produce a high- velocity isolated neutron 
star and an OB runaway star |38] . The high probability of disruption explains qualitatively why 
so few normal pulsars have companions. There are currently four known normal radio pulsars 
with massive main sequence companions in eccentric orbits which are examples of binary systems 
which survived the supernova explosion [1401 11531 12951 I^Sl I193| . Over the next 10 7-8 yr after the 
explosion, the neutron star may be observable as a normal radio pulsar spinning down to a period 
> several seconds. After this time, the energy output of the star diminishes to a point where it no 
longer produces significant amounts of radio emission. 

For those few binaries that remain bound, and in which the companion is sufficiently massive 
to evolve into a giant and overflow its Roche lobe, the old spun-down neutron star can gain a new 
lease of life as a pulsar by accreting matter and angular momentum at the expense of the orbital 
angular momentum of the binary system 0. The term "recycled pulsar" is often used to describe 
such objects. During this accretion phase, the X-rays produced by the liberation of gravitational 
energy of the infalling matter onto the neutron star mean that such a system is expected to be 
visible as an X-ray binary. Two classes of X-ray binaries relevant to binary and millisecond pulsars 
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Figure 7: Cartoon showing various evolutionary scenarios involving binary pulsars. 
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exist: neutron stars with high-mass or low-mass companions. Detailed reviews of the X-ray binary 
population, including systems likely to contain black holes, can be found elsewhere |33) . 




P b (days) 

Figure 8: Eccentricity versus orbital period for a sample of 21 low-mass binary pulsars which are not 
in globular clusters, with the triangles denoting three recently discovered systems \294\ - The solid 
line shows the median of the predicted relationship between orbital period and eccentricity ''2491. 
Dashed lines show 95% the confidence limit about this relationship. The dotted line shows Pb oc 
e 2 . Figure provided by Ingrid Stairs \294Jj using an adaptation of the orbital period- eccentricity 
relationship tabulated by Fernando Camilo. 

In a high-mass X-ray binary, the companion is massive enough that it might also explode as 
a supernova, producing a second neutron star. If the binary system is lucky enough to survive 
the explosion, the result is a double neutron star binary. At least five and probably eight such 
systems are now known, the original example being PSR B1913+16 |129) - a 59-ms radio pulsar 
which orbits its companion every 7.75 hr IH3ISI3]- In this formation scenario, PSR B1913+16 is 
an example of the older, first-born, neutron star that has subsequently accreted matter from its 
companion. 

For many years, no clear example was known where the second-born neutron star was observed 
as a pulsar. The discovery of the double pulsar J0737— 3039 |4*4l 1198) . where a 22.7-ms recycled 
pulsar "A" orbits a 2.77-s normal pulsar "B" every 2.4 hr, has now provided a dramatic confirma- 
tion of this evolutionary model in which we identify A and B as the first and second-born neutron 
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stars respectively. Just how many more observable double pulsar systems exist in our Galaxy is 
not clear. Although the population of double neutron star systems in general is reasonably well un- 
derstood (see Section l3.4.1[l . given that the lifetime of the second born pulsar is less than one tenth 
that of the recycled pulsar, and that its radio beam is likely to be much smaller (see Section l3.2.3[) . 
the prospects of ever finding more than a few 0737-like systems are rather low. 



I 1 1 1 1" T" 
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Figure 9: Companion mass versus orbital period for binary pulsars showing the whole sample 
where, in the absence of mass determinations, statistical arguments based on a random distribution 
of orbital inclination angles (see Section \4-4\ ) have been used to constrain the masses as shown 
(Panel a), and only those with well determined companion masses (Panel b). The dashed lines 
show the uncertainties in the predicted relation \3U7\j . This relationship indicates that as these 
systems finished a period of stable mass transfer due to Roche-lobe overflow, the size and hence 
period of the orbit was determined by the mass of the evolved secondary star. Figure provided by 
Marten van Kerkwijk [330]. 

The companion in a low-mass X-ray binary evolves and transfers matter onto the neutron star 
on a much longer time-scale, spinning it up to periods as short as a few ms |2]. Tidal forces during 
the accretion process serve to circularize the orbit. At the end of the spin-up phase, the secondary 
sheds its outer layers to become a white dwarf in orbit around a rapidly spinning millisecond 
pulsar. This model has gained strong support in recent years from the discoveries of quasi-periodic 
kHz oscillations in a number of low-mass X-ray binaries |338| . as well as Doppler-shifted 2.49-ms 
X-ray pulsations from the transient X-ray burster SAX J1808. 4-3658 pHOl 155] . Six other "X-ray 
millisecond pulsars" are now known with spin rates and orbital periods ranging between 185-600 Hz 
and 40 min-4.3 hr respectively [3391 12201 1144| . Despite intensive searches 03], n0 radio pulsations 
have so far been detected in these binaries. This could be a result of free-free absorption of any 
radio waves by the thick accretion disk, or perhaps quenching the accelerating potential in the 
neutron star magnetosphere by infalling matter. More sensitive radio observations are ultimately 
required to place stringent limits on the physical processes responsible for the lack of emission. 

Numerous examples of these systems in their post X-ray phase are now seen as the millisec- 
ond pulsar-white dwarf binary systems. Presently, 20 of these systems have compelling optical 
identifications of the white dwarf companion, and upper limits or tentative detections have been 
found in about 30 others |330| . Comparisons between the cooling ages of the white dwarfs and the 
millisecond pulsars confirm the age of these systems and suggest that the accretion rate during the 
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spin- up phase was well below the Eddington limit |114| . 

Further support for the above evolutionary scenarios comes from two correlations in the ob- 
served sample of low-mass binary pulsars. Firstly, as seen in Figure|Sl there is a strong correlation 
between orbital period and eccentricity. The data are very good agreement with a theoretical 
relationship which predicts a relic orbital eccentricity due to convective eddy currents in the ac- 
cretion process |249| . Secondly, as shown in Panel b of Figure El where companion masses have 
been measured accurately, through radio timing (see Section l4~4")l and/or through optical observa- 
tions |330| . they are in good agreement with a relation between companion mass and orbital period 
predicted by binary evolution theory |307| . A word of caution is required in using these models 
to make predictions, however. When confronted with a larger ensemble of binary pulsars using 
statistical arguments to constrain the companion masses (see Panel a of Figure EJ), current models 
have problems in explaining the full range of orbital periods on this diagram |294) . 

The range of white dwarf masses observed is becoming broader. Since this article originally 
appeared |179) . the number of "intermediate-mass binary pulsars" nas grown significantly |54) . 
These systems are distinct to the millisecond pulsar-white dwarf binaries in several ways: 

1. The spin period of the radio pulsar is generally longer (9-200 ms). 

2. The mass of the white dwarf is larger (typically > 0.5M Q ). 

3. The orbit, while still essentially circular, is often significantly more eccentric (e > 10~ 3 ). 

4. They do not necessarily follow the mass-period or eccentricity-period relationships. 

It is not presently clear whether these systems originated from low- or high-mass X-ray binaries. 
It was suggested by van den Heuvel EZQ] that they have more in common with high-mass systems. 
More recently, it has been proposed |172] that a thermal-viscous instability in the accretion disk of 
a low-mass X-ray binary could truncate the accretion phase and produce a more slowly spinning 
neutron star. 

2.7 Isolated recycled pulsars 

The scenarios outlined qualitatively above represent a reasonable understanding of binary evolu- 
tion. There are, however, a number of pulsars with spin properties that suggest a phase of recycling 
took place but have no orbiting companions. While the existence of such systems in globular clus- 
ters are more readily explained by the high probability of stellar interactions compared to the 
disk 283), it is somewhat surprising to find them in the Galactic disk. Out of a total of 66 millisec- 
ond pulsars in the Galactic disk, 15 are isolated (see Table|21). Although it has been proposed that 
these millisecond pulsars have ablated their companion via their strong relativistic winds |158j as 
may be happening in the PSR B1957+20 system |103| . it is not clear whether the energetics or 
time-scales for this process are feasible |170) . There is some observational evidence that suggests 
that solitary millisecond pulsars are less luminous than binary millisecond pulsars |23l I165| . If 
confirmed by future discoveries, this would need to be explained by any viable evolutionary model. 

There are two further "anomalous" isolated pulsars with periods in the range 55-60 ms [5511188) . 
When placed on the P-P diagram, these objects populate the region occupied by the double 
neutron star binaries. The most natural explanation for their existence, therefore, is that they 
are "failed double neutron star binaries" which disrupted during the supernova explosion of the 
secondary |SSJ. A simple calculation |188) . however, suggests that for every double neutron star 
we should see of order ten such isolated objects. Although some other examples of such pulsars 
are known |150l I9l|l42| . exactly why so few are observed is currently not clear. 
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2.8 Pulsar velocities 



Pulsars have long been known to have space velocities at least an order of magnitude larger than 
those of their main sequence progenitors, which have typical values between 10 and 50 km s _1 . 
The first direct evidence for large velocities came from optical observations of the Crab pulsar [32 1) . 
showing that the neutron star has a velocity in excess of 100 km s _1 . Proper motions for 233 pulsars 
have subsequently been measured largely by radio timing and interferometric techniques [1941 1241 
I9b1lll5lll22) . These data imply a broad velocity spectrum ranging from to over 1000 km s" 1 [2UT] . 



Figure 10: GIF movie showing a simulation following the motion of 100 pulsars in a model gravita- 
tional potential of our Galaxy for 200 Myr. The view is edge-on, i.e. the horizontal axis represents 
the Galactic plane (30 kpc across) while the vertical axis represents ±10 kpc from the plane. This 
snapshot shows the initial configuration of young neutron stars. 

Such large velocities are perhaps not surprising, given the violent conditions under which neu- 
tron stars are formed. Shklovskii |282) demonstrated that, if the explosion is only slightly asym- 
metric, an impulsive "kick" velocity of up to 1000 km s _1 can be imparted to the neutron star. In 
addition, if the neutron star progenitor was a member of a binary system prior to the explosion, 
the pre-supernova orbital velocity will also contribute to the resulting speed of the newly-formed 
pulsar. As Figure QJU illustrates, high-velocity pulsars born close to the Galactic plane quickly 
migrate to higher Galactic latitudes. Given such a broad velocity spectrum, as many as half of all 
pulsars will eventually escape the Galactic gravitational potential [2011 171| . 

The distribution of pulsar velocities remains an issue of contention. For example, Monte Carlo 
simulations by Arzoumanian et al. 8 strongly favour a bimodal distribution with low and high 
velocity peaks. On the other hand, a recent study (122) found the mean birth velocity of normal 
pulsars to be consistent with a Maxwellian distribution with a mean of ~ 400 km s _1 . Regardless 
of the form of the distribution, however, we can say that the mean velocity of young pulsars is 
significantly larger than the 80-140 km s _1 range for millisecond and binary pulsars [1781 1701 |205j . 
The main reasons for the lower velocities are the fact that the kick must have been small to avoid 
disruption, and the surviving neutron star has to pull the companion along with it, thus slowing 
the system down. 
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2.9 Current and future binary and millisecond pulsar search strategies 

2.9.1 All-sky searches 

The oldest radio pulsars form a relaxed population of stars oscillating in the Galactic gravitational 
potential [113) . The scale height for such a population is at least 500 pc, about 10 times that of 
the massive stars which populate the Galactic plane. Since the typical ages of millisecond pulsars 
are several Gyr or more, we expect, from our vantage point in the Galaxy, to be in the middle of an 
essentially isotropic population of nearby sources. All- sky searches for millisecond pulsars at high 
Galactic latitudes have been very effective in probing this population. Much of the initial interest 
and excitement in this area began with the discovery of two recycled pulsars at high latitudes with 
the 305-m Arecibo telescope: the double neutron star binary B1534+12 |343j and the "planets 
pulsar" B1257+12 |346) . Surveys carried out at Arecibo, Parkes, Jodrell Bank and Green Bank 
(using the 140 ft telescope) by others in the 1990s (summarised in three review papers |4*71 ISTfl 1ST] ) 
have found many other millisecond and recycled pulsars in this way. 

2.9.2 Searches close to the plane of our Galaxy 

Young pulsars are most likely to be found near to their places of birth, before they have had time 
for their velocity to move them away, and hence they lie close to the Galactic plane (see Figure EJ. 
This was the target region of the main Parkes multibeam survey and has already resulted in 
the discovery of over 700 new pulsars |210l 12211 11631 11231 194) . almost half the number currently 
known! Such a large haul inevitably results in a number of interesting individual objects such 
as PSR J1141— 6545, a young pulsar in a relativistic 4.8-hr orbit around a white dwarf [1541 124 ll 
I2"51ll27) . PSR J1740— 3052, a young pulsar orbiting an ~ 11 M@ star (probably a main sequence 
B-star [2951 1296) ). PSR 1638—4715, a young pulsar in a ~ 5 yr-eccentric orbit (e ~ 0.9) around 
around a 10-20 M companion |193| . several intermediate-mass binary pulsars |S3], and two double 
neutron star binaries |199l I93| . Although the main survey has now been completed, extensions of 
the survey region, and re-analyses of existing data [HI] will ensure further discoveries in the near 
future. 

2.9.3 Searches at intermediate Galactic latitudes 

To probe more deeply into the population of millisecond and recycled pulsars than possible at 
high Galactic latitudes, the Parkes multibeam system was also used to survey intermediate lati- 
tudes H22ED1- Among the 69 new pulsars found in the survey, 8 are relatively distant recycled 
objects. Two of the new recycled pulsars from this survey |90| are mildly relativistic neutron star- 
white dwarf binaries. An analysis of the full results from this survey should significantly improve 
our knowledge on the Galaxy-wide population and birth-rate of millisecond pulsars. Arecibo sur- 
veys at intermediate latitudes also continue to find new pulsars, such as the long-period binaries 
J2016+1948 and J0407+1607 [2"3T1 IT977] . and the likely double neutron star system J1829+2456 (H2J 

2.9.4 Targeted searches of globular clusters 

Globular clusters have long been known to be breeding grounds for millisecond and binary pul- 
sars The main reason for this is the high stellar density and consequently high rate of stellar 
interaction in globular clusters relative to most of the rest of the Galaxy. As a result, low-mass 
X-ray binaries are almost 10 times more abundant in clusters than in the Galactic disk. In addi- 
tion, exchange interactions between binary and multiple systems in the cluster can result in the 
formation of exotic binary systems [283 . To date, searches have revealed 103 pulsars in 24 globular 
clusters (see Table [5] and 230 58 ) . Early highlights include the double neutron star binary in 
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M15 |257) and a low-mass binary system with a 95-min orbital period in 47 Tucanae , one of 
22 millisecond pulsars currently known in this cluster alone [S3] H83) . 

On-going surveys of clusters continue to yield new surprises [2641 177) including the discovery 
of the most eccentric binary pulsar so far - J0514— 4002 is a 4.99 ms pulsar in a highly eccentric 
(e = 0.89) binary system in the globular cluster NGC 1851 jlOO) . The latest sensation, however, is 
the discovery |266l 1233] of 28 millisecond pulsars in Terzan 5 with the Green Bank Telescope |232) . 
This brings the total known in this cluster to 26. The spin periods and orbital parameters of the 
new pulsars reveal that, as a population, they are significantly different to the pulsars of 47 Tucanae 
which have periods in the range 2-8 ms |183j . The spin periods of the new pulsars span a much 
broader range (1.67-80 ms) including the second and third shortest spin periods of all pulsars 
currently known. The binary pulsars include two systems with eccentric orbits and likely white 
dwarf companions. No such systems are known in 47 Tucanae. The difference between the two 
pulsar populations may reflect the different evolutionary states and physical conditions of the two 
clusters. In particular, the central stellar density of Terzan 5 is about twice that of 47 Tucanae, 
suggesting that the increased rate of stellar interactions might disrupt the recycling process for the 
neutron stars in some binary systems and induce larger eccentricities in others. 

2.9.5 Current and future surveys 

Motivated by the successes at Parkes, a multibeam receiver has recently been installed at the 
Arecibo telescope |227) . A preliminary survey using this instrument began in 2004 and has so 
far discovered 10 pulsars 229 . The main virtue of this survey is the shorter integration times 
employed and resulting higher sensitivity to accelerated systems than the Parkes survey. Depending 
on the stamina of the observers, up to 1000 pulsars could be found in this survey over the next 
decade. Other surveys with the Green Bank Telescope |232| and with the Giant Metre Wave Radio 
Telescope |231) are also expected to make major contributions. 

All surveys that have so far been conducted, or will be carried out in the next few years, are 
likely to be surpassed by the Square Kilometre Array |132| - the next generation radio telescope, 
planned to come online around 2020. Simulations suggest j!62| that at least 10 4 pulsars, including 
10 3 millisecond pulsars, could be detected in our Galaxy. 

2.10 Going further 

In this brief review of pulsar astronomy, we have only skimmed the surface of many topics which 
are covered more thoroughly elsewhere. Two recent textbooks, Pulsar Astronomy 200 and Hand- 
book of Pulsar Astronomy 185 , cover the literature and techniques and provide excellent further 
reading. The morphological properties of pulsars have recently been comprehensively discussed 
in two recent reviews 12861 1280) . Those wishing to approach the subject from a more theoret- 
ical viewpoint are advised to read The Theory of Neutron Star Magnetospheres |218) and The 
Physics of the Pulsar Magnetosphere |32| . Our summary of evolutionary aspects serves merely as 
a primer to the vast body of literature available. Further insights can be found from more detailed 
reviews [331 12311 1232] . 

Pulsar resources available on the Internet are continually becoming more extensive and useful. 
A good starting point for pulsar-surfers is the Handbook of Pulsar Astronomy website |186) . as 
well as the pages maintained at Arecibo |226) . Berkeley |325) . Bonn |215| . Cagliari |324) . Jodrcll 
Bank 132^], Princeton [233], Swinburne [305], UBC [323] and Sydney [255] . 
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3 Pulsar Statistics and Demography 



The observed pulsar sample is heavily biased towards the brighter objects that are the easiest 
to detect. What we observe represents only the tip of the iceberg of a much larger underlying 
population |108| . The extent to which the sample is incomplete is well demonstrated by the 
projection of pulsars onto the Galactic plane and their cumulative number distribution as a function 
of distance shown in Figure 1111 The clustering of sources around the Sun seen in the left panel 
of Figure 1111 is clearly at variance with the distribution of other stellar populations which show a 
radial distribution about the Galactic centre. 

The extent to which the pulsar sample is incomplete can be seen from the cumulative number of 
pulsars as a function of the projected distance from the Sun. In FigurelTDthe observed distribution 
is compared to the expected distribution for a simple model population with no selection effects. 
The observed number distribution becomes strongly deficient beyond a few kpc. 




-10 10 0.01 0.1 1 10 
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Figure 11: Left panel: The current sample of all known radio pulsars projected onto the Galactic 
plane. The Galactic centre is at the origin and the Sun is at (0,8.5) kpc. Note the spiral-arm 
structure seen in the distribution which is now required by the electron density model \7S\ \74\j - 
Right panel: Cumulative number of pulsars as a function of projected distance from the Sun. The 
solid line shows the observed sample while the dotted line shows a model population free from 
selection effects. 



3.1 Selection effects in pulsar searches 
3.1.1 The inverse square law and survey thresholds 

The most prominent selection effect is the inverse square law, i.e. for a given intrinsic luminosity 1 , 
the observed flux density falls off as the inverse square of the distance. This results in the observed 
sample being dominated by nearby and/or high luminosity objects. Beyond distances of a few kpc 

1 Pulsar astronomers usually define the luminosity L = Sd? , where S is the mean flux density at 400 MHz (a 
standard observing frequency) and d is the distance derived from the DM (see Section 12.41 . Since this ignores any 
assumptions about beaming or geometrical factors, it is sometimes referred to as a "pseudoluminosity" 0. 



15 



from the Sun, the apparent flux density falls below the detection thresholds S m in of most surveys. 
Following J53J, we express this threshold as follows: 

- - m m wr m~" 2 (¥r w" • 

where S/N m ; n is the threshold signal-to-noise ratio, 77 is a generic fudge factor (< 1) which ac- 
counts for losses in sensitivity (e.g., due to sampling and digitization noise), n po i is the number of 
polarizations recorded (either 1 or 2), T roc and T s k y are the receiver and sky noise temperatures, 
G is the gain of the antenna, Av is the observing bandwidth, ij nt is the integration time, W is the 
detected pulse width and P is the pulse period. 



3.1.2 Interstellar pulse dispersion and multipath scattering 

It follows from Equation © that the sensitivity decreases as W/(P — W) and hence W increases. 
Also note that if W > P, the pulsed signal is smeared into the background emission and is no 
longer detectable, regardless of how luminous the source may be. The detected pulse width W may 
be broader than the intrinsic value largely as a result of pulse dispersion and multipath scattering 
by free electrons in the interstellar medium. The dispersive smearing scales as Ai//i/ 3 , where v is 
the observing frequency. This can largely be removed by dividing the pass-band into a number of 
channels and applying successively longer time delays to higher frequency channels before summing 
over all channels to produce a sharp profile. This process is known as incoherent dedispersion. 

The smearing across the individual frequency channels, however, still remains and becomes 
significant at high dispersions when searching for short-period pulsars. Multipath scattering from 
electron density irregularities results in a one-sided broadening due to the delay in arrival times. 
A simple scattering model is shown in Figure [T2l in which the scattering electrons are assumed to 
lie in a thin screen between the pulsar and the observer 278 . The timescale of this effect varies 
roughly as i^ -4 , which can not be removed by instrumental means. 
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Figure 12: Left panel: Pulse scattering caused by irregularities in the interstellar medium. The 
different path lengths and travel times of the scattered rays result in a "scattering tail" in the 
observed pulse profile which lowers its signal-to-noise ratio. Right panel: A simulation showing 
the percentage of Galactic pulsars that are likely to be undetectable due to scattering as a function 
of observing frequency. Low-frequency (< 1 GHz) surveys clearly miss a large percentage of the 
population due to this effect. 

Dispersion and scattering are most severe for distant pulsars in the inner Galaxy where the 
number of free electrons along the line of sight becomes large. The strong frequency dependence of 
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both effects means that they are considerably less of a problem for surveys at observing frequencies 
> 1400 MHz |67II139| compared to the 400-MHz search frequency used in early surveys. An added 
bonus for such observations is the reduction in T s k y , since the spectral index of the non-thermal 
Galactic emission is about —2.8 |169| . Pulsars themselves have typical spectral indices of — 1.6 |191) . 
so that flux densities are roughly an order of magnitude lower at 1400 MHz compared to 400 MHz. 
Fortunately this can be at least partially compensated for by the use of larger receiver bandwidths 
at higher radio frequencies. For example, the 1370-MHz system at Parkes has a bandwidth of 
288 MHz [T39"] compared to the 430-MHz system, where nominally 32 MHz is available [3TT] . 

3.1.3 Orbital acceleration 

Standard pulsar searches use Fourier techniques |185| to search for a-priori unknown periodic 
signals and usually assume that the apparent pulse period remains constant throughout the ob- 
servation. For searches with integration times much greater than a few minutes, this assumption 
is only valid for solitary pulsars or binary systems with orbital periods longer than about a day. 
For shorter-period binary systems, the Doppler-shifting of the period results in a spreading of the 
signal power over a number of frequency bins in the Fourier domain, leading to a reduction in 
S/N |13fij . An observer will perceive the frequency of a pulsar to shift by an amount aT/(Pc), 
where a is the (assumed constant) line-of-sight acceleration during the observation of length T, P 
is the (constant) pulsar period in its rest frame and c is the speed of light. Given that the width 
of a frequency bin in the Fourier domain is 1/T, we see that the signal will drift into more than 
one spectral bin if aT 2 /(Pc) > 1. Survey sensitivities to rapidly-spinning pulsars in tight orbits 
are therefore significantly compromised when the integration times are large. 

As an example of this effect, as seen in the time domain, Figure ITTA shows a 22.5-min search 
mode observation of the binary pulsar B1913+16 |129ll312ll313) . Although this observation covers 
only about 5% of the orbit (7.75 hr), the severe effects of the Doppler smearing on the pulse signal 
are very apparent. While the standard search code nominally detects the pulsar with S/N = 9.5 
for this observation, it is clear that this value is significantly reduced due to the Doppler shifting 
of the pulse period seen in the individual sub-integrations. 

It is clearly desirable to employ a technique to recover the loss in sensitivity due to Doppler 
smearing. One such technique, the so-called "acceleration search" |219) . assumes the pulsar has 
a constant acceleration during the observation. Each time series can then be re-sampled to refer 
it to the frame of an inertial observer using the Doppler formula to relate a time interval r in 
the pulsar frame to that in the observed frame at time t, as r(t) oc (1 + at/c). Searching over 
a range of accelerations is desirable to find the time series for which the trial acceleration most 
closely matches the true value. In the ideal case, a time series is produced with a signal of constant 
period for which full sensitivity is recovered (see right panel of Figure 1131) . This technique was 
first used to find PSR B2127+11C 0, a double neutron star binary in M15 which has parameters 
similar to B1913+16. More recently, its application to 47 Tucanae jSS] resulted in the discovery 
of nine binary millisecond pulsars, including one in a 96-min orbit around a low-mass (0.15 M©) 
companion. This is currently the shortest binary period for any known radio pulsar. 

For intermediate orbital periods, in the range 30 min-several hours, another promising tech- 
nique is the dynamic power spectrum search. Here the time series is split into a number of smaller 
contiguous segments which are Fourier-transformed separately. The individual spectra are dis- 
played as a two-dimensional (frequency versus time) image. Orbitally modulated pulsar signals 
appear as sinusoidal signals in this plane as shown in Figure ITU 

This technique has been used by various groups where spectra are inspected visually |206) . 
Much of the human intervention can be removed using a hierarchical scheme for selecting significant 
events jHS]- This approach was recently applied to a search of the globular cluster M62 resulting 
in the discovery of three new pulsars. One of the new discoveries - M62F, a faint 2.3-ms pulsar in 
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Figure 13: Left panel: A 22.5-min Arecibo observation of the binary pulsar B1913+16. The as- 
sumption that the pulsar has a constant period during this time is clearly inappropriate given the 
quadratic drifting in phase of the pulse during the observation (linear grey scale plot). Right panel: 
The same observation after applying an acceleration search. This shows the effective recovery of 
the pulse shape and a significant improvement in the signal-to-noise ratio. 
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a 4.8-hr orbit - was detectable only using the dynamic power spectrum technique. 

For the shortest orbital periods, the assumption of a constant acceleration during the ob- 
servation clearly breaks down. In this case, a particularly efficient algorithm has been devel- 
oped |262l 11431 1263| which is optimised to finding binaries with periods so short that many orbits 
can take place during an observation. This "phase modulation" technique exploits the fact that 
the Fourier components are modulated by the orbit to create a family of periodic sidebands around 
the nominal spin frequency of the pulsar. While this technique has so far not resulted in any new 
discoveries, the existence of short period binaries in 47 Tucanae Terzan 5 | 266| and the 11-min 
X-ray binary X1820— 303 in NGC 6624 |301) . suggests that there may be more ultra-compact binary 
pulsars that await discovery. 

3.2 Correcting the observed pulsar sample 

In the following, we review common techniques to account for the various selection effects and form 
a less biased picture of the true pulsar population. 

3.2.1 Scale factor determination 

A very useful technique [2501 1334] . is to define a scaling factor £ as the ratio of the total Galactic 
volume weighted by pulsar density to the volume in which a pulsar is detectable: 



Here, £(-R, z) is the assumed pulsar space density distribution in terms of galactocentric radius 
R and height above the Galactic plane z. Note that £ is primarily a function of period P and 
luminosity L such that short-period/low-luminosity pulsars have smaller detectable volumes and 
therefore higher £ values than their long-period/high-luminosity counterparts. 

In practice, £ is calculated for each pulsar separately using a Monte Carlo simulation to model 
the volume of the Galaxy probed by the major surveys |222| . For a sample of iV b s observed pulsars 
above a minimum luminosity £ m i n , the total number of pulsars in the Galaxy is 



where / is the model-dependent "beaming fraction" discussed below in Section 2. 31 Note that 
this estimate applies to those pulsars with luminosities > i m i n - Monte Carlo simulations have 
shown this method to be reliable, as long as -/V b s is reasonably large |182j . 

3.2.2 The small-number bias 

For small samples of observationally-selected objects, the detected sources are likely to be those 
with larger-than-average luminosities. The sum of the scale factors (JSJ), therefore, will tend to 
underestimate the true size of the population. This "small-number bias" was first pointed out [1451 
1151) for the sample of double neutron star binaries where we know of only four systems relevant 
for calculations of the merging rate (see Section l3.4.1fl . Only when iV bs > 10 does the sum of the 
scale factors become a good indicator of the true population size. 

Despite a limited sample size, recent work |156) has demonstrated that rigorous confidence 
intervals of Nq can be derived. Monte Carlo simulations verify that the simulated number of 
detected objects A^etected closely follows a Poisson distribution and that iVdetected = o>Nq, where a 
is a constant. By varying the value of Nq in the simulations, the mean of this Poisson distribution 



J J Q Y,(R,z)RdRdz 
J J PL E{R, z)RdRdz' 



(4) 




(5) 



19 



Model Galactic Population 




Figure 15: Small-number bias of the scale factor estimates derived from a synthetic population 
of sources where the true number of sources is known. Left panel: An edge-on view of a model 
Galactic source population. Right panel: The thick line shows Nq, the true number of objects in 
the model Galaxy, plotted against iV bservedj the number detected by a flux-limited survey. The 
thin solid line shows N est , the median sum of the scale factors, as a function of N ^ s from a large 
number of Monte Carlo trials. Dashed lines show 25 and 75% percentiles of the N est distribution. 

can be measured. Using a Bayesian analysis it was shown |156| that, for a single object, the 
probability density function of the total population is 

P(N G ) = a 2 Nq exp(-a7V G )- (6) 

Adopting the necessary assumptions required in the Monte Carlo population about the underlying 
pulsar distribution, this technique can be used to place interesting constraints on the size and, as 
we shall see later, birth rate of the underlying population. 

3.2.3 The beaming correction 

The "beaming fraction" / in Equation (JSJ is the fraction of Air steradians swept out by a pulsar's 
radio beam during one rotation. Thus / is the probability that the beam cuts the line-of-sight 
of an arbitrarily positioned observer. A naive estimate for / of roughly 20% assumes a circular 
beam of width ^10° and a randomly distributed inclination angle between the spin and magnetic 
axes |311) . Observational evidence suggests that shorter period pulsars have wider beams and 
therefore larger beaming fractions than their long-period counterparts |223l 12021 1341 13(H)] . As can 
be seen in Figure however, a consensus on the beaming fraction-period relation has yet to be 
reached. 

When most of these beaming models were originally proposed, the sample of millisecond pul- 
sars was < 5 and hence their predictions about the beaming fractions of short-period pulsars 
relied largely on extrapolations from the normal pulsars. An analysis of a large sample of mil- 
lisecond pulsar profiles |165) suggests that their beaming fraction lies between 50 and 100%. The 
large beaming fraction and narrow pulses often observed strongly suggests a fan beam model for 
millisecond pulsars |218) . 
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Figure 16: Beaming fraction plotted against pulse period for four different beaming models: Narayan 
& Vivekanand 1983 (NV83) JMUtl , Lyne & Manchester 1988 (LM88) '2021, Biggs 1990 (JDB90) 
and Tauris & Manchester 1998 (TM88) TMb^ . 



3.3 The population of normal and millisecond pulsars 
3.3.1 Luminosity distributions and local number estimates 

Based on a number of all-sky surveys carried out in the 1990s, the scale factor approach has been 
used to derive the characteristics of the true normal and millisecond pulsar populations and is 
based on the sample of pulsars within 1.5 kpc of the Sun |2()5) . Within this region, the selection 
effects are well understood and easier to quantify than in the rest of the Galaxy. These calculations 
should therefore give a reliable local pulsar population estimate. 

The luminosity distributions obtained from this analysis are shown in Figure 1171 For the 
normal pulsars, integrating the corrected distribution above 1 mJy kpc 2 and dividing by n x 
(1.5) 2 kpc 2 yields a local surface density, assuming a beaming model of 156± 31pulsars kpc~ 2 
for luminosities above 1 mJy kpc 2 . The same analysis for the millisecond pulsars, assuming a mean 
beaming fraction of 75% |165) . leads to a local surface density of 38 ± 16 pulsars kpc -2 also for 
luminosities above 1 mJy kpc 2 . 



3.3.2 Galactic population and birth-rates 

Integrating the local surface densities of pulsars over the whole Galaxy requires a knowledge of the 
presently rather uncertain Galactocentric radial distribution [13711181] . One approach is to assume 
that pulsars have a radial distribution similar to that of other stellar populations and to scale the 
local number density with this distribution in order to estimate the total Galactic population. 
The corresponding local-to-Galactic scaling is 1000 ± 250 kpc 2 |268j . This implies a population of 
~ 160, 000 active normal pulsars and ~ 40, 000 millisecond pulsars in the Galaxy. 

Based on these estimates, we are in a position to deduce the corresponding rate of formation 
or birth-rate required to sustain the observed population. From the P-P diagram in Figure [3] 
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Figure 17: Left panel: The corrected luminosity distribution (solid histogram with error bars) for 
normal pulsars. The corrected distribution before the beaming model has been applied is shown by 
the dot-dashed line. Right panel: The corresponding distribution for millisecond pulsars. In both 
cases, the observed distribution is shown by the dashed line and the thick solid line is a power law 
with a slope of — 1. The difference between the observed and corrected distributions highlights the 
severe under-sampling of low-luminosity pulsars. 

we infer a typical lifetime for normal pulsars of ~ 10 7 yr, corresponding to a Galactic birth rate 
of ~ 1 per 60 yr - consistent with the rate of supernovae jH2Hj- As noted in Section l2~2l the 
millisecond pulsars are much older, with ages close to that of the Universe r u (we assume here 
t u = 13.8 Gyr [347] ) . Taking the maximum age of the millisecond pulsars to be t u , we infer a 
mean birth rate of at least 1 per 345, 000 yr. This is consistent, within the uncertainties, with the 
birth-rate of low-mass X-ray binaries |189) . 

3.4 The population of relativistic binaries 

Although no radio pulsar has so far been observed in orbit around a black hole companion, we 
now know of several double neutron star and neutron star-white dwarf binaries which will merge 
due to gravitational wave emission within a reasonable timescale. The current sample of objects 
is shown as a function of orbital period and eccentricity in Figure IT%1 Isochrones showing various 
coalescence times r g are calculated using the expression 

where mi and mi are the masses of the two stars, [i = m\m%j{rn\ + mi) is the so-called "reduced 
mass" , Pb is the binary period and e is the eccentricity. This formula is a good analytic approxi- 
mation (within a few percent) to the numerical solution of the exact equations for r g 2451 1246] . 

In addition to testing general relativity through observations of these systems (see Section |4~4"|) . 
estimates of their Galactic population and merger rate are of great interest as one of the prime 
sources for current gravitational wave detectors such as GEO600 |112) . LIGO $jg\, VIRGO |131) 
and TAMA p25 . In the following, we review recent empirical determination of the population 
sizes and merging rates of binaries where at least one component is visible as a radio pulsar. 
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Figure 18: The relativistic binary merging plane. Top: Orbital eccentricity versus period for eccen- 
tric binary systems involving neutron stars. Bottom: Orbital period distribution for the massive 
white dwarf-pulsar binaries. 



3.4.1 Double neutron star binaries 

As discussed in Section |2~21 double neutron star (DNS) binaries are expected to be rare. This is 
certainly the case; despite extensive searches, only five certain DNS binaries are currently known: 
PSRs J0737-3039 g2], B1534+12 [313], J1756-2251 [H2J, B1913+16 [T2J| and B2127+11C [21771 ■ 
Although we only see both neutron stars as pulsars in J0737— 3039 |198) . we are "certain" of 
the identification in the other four systems from precise component mass measurements from 
pulsar timing observations (see Section IPjl . The spin and orbital parameters, timescales and 
mass constraints for these systems are listed in Table ^ along with three other binaries with 
eccentric orbits, mass functions and periastron advance measurements that are consistent with a 
DNS identification, but for which there is presently not sufficient component mass information to 
confirm their nature. 

Despite the uncertainties in identifying DNS binaries, for the purposes of determining the 
Galactic merger rate, the systems 2 for which r g is less than r u (i.e. PSRs J0737— 3039, B1534+12, 
J1756-2251, B1913+16 and B2127+11C) are primarily of interest. Of these PSR B2127+11C is 
in the process of being ejected from the globular cluster M15 [2571 1252) and is thought to make 
only a negligible contribution to the merger rate |248) . The general approach with the remaining 
systems is to derive scale factors for each object, construct the probability density function of 
their total population (as outlined in Section 13. 2.1 [I and then divide these by a reasonable estimate 
for the lifetime. It is generally accepted |151j that the observable lifetimes for these systems are 
determined by the timescale on which the current orbital period is reduced by a factor of two jjJJ . 
Below this point, the orbital smearing selection effect discussed in Section f3 . 1 . 31 will render the 

2 The likely DNS binary J1756— 2251 1931 has a merger time of 1.7 Gyr and will make a small contribution to the 
rate. It has so far not been incorporated into the latest calculations. 



23 





J0737-3039 


J1518+4904 


B1534+12 


J1756- 


2251 


/ > r J] 

P [msj 


22.1 / 277\) 


40.9 


o t a 

37.9 




z8.5 


ft [dj 


r\ 1 no 


q a. 


U.4 




U.oz 


e 


A AO O 

0.088 


0.25 


A 

0.27 




A 1 O 

0.18 


r c [10° yrj 


210/50 


200 


2.5 




A A 

4.4 




Q7 
O ( 


z4UUU 


Z 1 




1D9 


Masses measured? 


Yes 


No 


Yes 




Yes 




J1811-1736 


J1829+2456 


B1913+16 


B2127- 


MIC 


P [ms] 


104.2 


41.0 


59.0 




30.5 


ft [d] 


18.8 


1.18 


0.3 




0.3 


e 


0.83 


0.14 


0.62 




0.68 


Tc [10 8 yr] 


9.2 


130 


1.1 




0.97 


r g [10 8 yr] 


10000 


600 


3.0 




2.2 


Masses measured? 


No 


No 


Yes 




Yes 



Table 1: Known and likely DNS binaries. Listed are the pulse period P , orbital period Pb, orbital 
eccentricity e, characteristic age r c , and expected binary coalescence time-scale r g due to gravita- 
tional wave emission calculated from Equation Qj. To distinguish between definite and candidate 
DNS systems, we also list whether the masses of both components have been determined via the 
measurement of two or more post Keplerian parameters as described in Section\4-4\ 



binary undetectable by current surveys. 

The results of the most recent study of this kind [14711148] . which take into account the discovery 
of PSR J0737— 3039, are summarised in the left panel of Figure H^l The combined Galactic merger 
rate, dominated by the double pulsar, is found to be 80^70° Myr -1 , where the uncertainties reflect 
the 95% confidence level using the techniques summarised in Section 13.2.21 Extrapolating this 
number to include DNS binaries detectable by LIGO in other galaxies |248) . the expected event 
rate is 35lgg x 10~ 3 yr -1 for initial LIGO and 190^150 y r_1 f° r advanced LIGO. Future prospects 
for detecting gravitational wave emission from binary neutron star inspirals are therefore very 
encouraging. Since much of the uncertainty in the rate estimates is due to our ignorance of the 
underlying distribution of double neutron star systems, future detection rates could ultimately 
constrain the properties of this exciting binary species. 

Although the double pulsar system J0737— 3039 will not be important for ground-based de- 
tectors until its final coalescence in another 85 Myr, it may be a useful calibration source for the 
future space-based detector LISA |224| . Recent calculations |146| show that a 1-yr observation 
with LISA would detect (albeit with S/N ~ 2) the continuous emission at a frequency of 0.2 mHz 
based on the current orbital parameters. Although there is the prospect of using LISA to detect 
similar systems systems through their continuous emission, current calculations |146) suggest that 
significant (S/N > 5) detections are not likely. Despite these limitations, it is likely that LISA 
observations will be able to place independent constraints on the Galactic DNS binary population 
after several years of operation. 

3.4.2 White dwarf neutron star binaries 

Although the population of white dwarf-neutron star (WDNS) binaries in general is substantial, 
the fraction which will merge due to gravitational wave emission is small. Like the DNS binaries, 
the observed WDNS sample suffers from small-number statistics. From Figure ^0 we note that 
only three WDNS systems are currently known that will merge within r u , PSRs J0751+1807 192 , 
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Figure 19: The current best empirical estimates of the coalescence rates of relativistic binaries 
involving neutron stars. The individual contributions from each known binary system are shown as 
dashed lines, while the solid line shows the total probability density function on a logarithmic and 
(inset) linear scale. The left panel shows DNS binaries \H1\ while the right panel shows the 

equivalent results for NS~ WD binaries \157\[ . 
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J1757— 5322 "OH] an d J1141— 6545 |154) . Applying the same techniques as used for the DNS popu- 
lation, the merging rate contributions of the three systems can be calculated |157j and are shown 
in Figure ED The combined Galactic coalescence rate is 4+| Myr -1 (at 68% confidence interval). 
Although the orbital frequencies of these objects at coalescence are too low to be detected by 
LIGO, they do fall within the band planned for LISA |224) . Unfortunately, an extrapolation of 
the Galactic event rate out to distances at which such events would be detectable by LISA does 
not suggest that these systems will be a major source of detection |157| . Similar conclusions were 
reached by considering the statistics of low-mass X-ray binaries [6*9] . 

3.5 Going further 

Studies of pulsar population statistics represent a large proportion of the pulsar literature. Good 
starting points for further reading can be found in other review articles [2181 13.3) . Our coverage of 
compact object coalescence rates has concentrated on empirical methods. An alternative approach 
is to undertake a full-blown Monte Carlo simulation of the most likely evolutionary scenarios 
described in Section l2.6l In this approach, a population of primordial binaries is synthesized with 
a number of underlying distribution functions: primary mass, binary mass ratio, orbital period 
distribution etc. The evolution of both stars is then followed to give a predicted sample of binary 
systems of all the various types. Although selection effects are not always taken into account in 
this approach, the final census is usually normalized to the star formation rate. 

Numerous population syntheses (most often to populations of binaries where one or both mem- 
bers are NSs) can be found in the literature [861 12721 13221 12541 11731 127) . A group in Moscow has 
made a web interface to their code 302 . Although extremely instructive, the uncertain assump- 
tions about initial conditions, the physics of mass transfer and the kicks applied to the compact 
object at birth result in a wide range of predicted event rates which are currently broader than the 
empirical methods [1481 11571 1149) . Ultimately, the detection statistics from the gravitational wave 
detectors could provide far tighter constraints on the DNS merging rate than the pulsar surveys 
from which these predictions are made. 
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4 Principles and Applications of Pulsar Timing 



It became clear soon after their discovery that pulsars are excellent celestial clocks. The period 
of the first pulsar |II9| was found to be stable to one part in 10 7 over a few months. Following 
the discovery of the millisecond pulsar B1937+21 ^J] it was demonstrated that its period could 
be measured to one part in 10 13 or better |84| . This unrivaled stability leads to a host of applica- 
tions including uses as time keepers, probes of relativistic gravity and natural gravitational wave 
detectors. 



4.1 Observing basics 

Each pulsar is typically observed at least once or twice per month over the course of a year 
to establish its basic properties. Figure |2"U1 summarises the essential steps involved in a "time-of- 
arrival" (TOA) measurement. Pulses from the neutron star traverse the interstellar medium before 
being received at the radio telescope where they are de-dispersed and added to form a mean pulse 
profile. 
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Figure 20: Schematic showing the main stages involved in pulsar timing observations. 

During the observation, the data regularly receive a time stamp, usually based on a caesium 
time standard or hydrogen maser at the observatory plus a signal from the Global Positioning 
System of satellites (GPS; see [EB])- The TOA is defined as the arrival time of some fiducial point 
on the integrated profile with respect to either the start or the midpoint of the observation. Since 
the profile has a stable form at any given observing frequency (see Section f2.3f) . the TOA can 
be accurately determined by cross-correlation of the observed profile with a high S/N "template" 
profile obtained from the addition of many observations at the particular observing frequency. 

Successful pulsar timing requires optimal TOA precision which largely depends on the signal- 
to-noise ratio (S/N) of the pulse profile. Since the TOA uncertainty exoA is roughly the pulse 
width divided by the S/N, using Equation @ we may write the fractional error as 

Here, S paT is the flux density of the pulsar, T rcc and T s ^y are the receiver and sky noise temperatures, 
G is the antenna gain, Av is the observing bandwidth, t int is the integration time, W is the pulse 
width and P is the pulse period (we assume W < P). Optimal results are thus obtained for 
observations of short period pulsars with large flux densities and small duty cycles (W/P) using 
large telescopes with low-noise receivers and large observing bandwidths. 
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One of the main problems of employing large bandwidths is pulse dispersion. As discussed 
in Section l2~4l pulses emitted at lower radio frequencies travel slower and arrive later than those 
emitted at higher frequencies. This process has the effect of "stretching" the pulse across a finite 
receiver bandwidth, increasing W and therefore increasing exoA- For normal pulsars, dispersion 
can largely be compensated for by the incoherent de-dispersion process outlined in Section l3~Tl 

The short periods of millisecond pulsars offer the ultimate in timing precision. In order to fully 
exploit this, a better method of dispersion removal is required. Technical difficulties in building 
devices with very narrow channel bandwidths require another dispersion removal technique. In the 
process of coherent de-dispersion [1101 1185) the incoming signals are de-dispersed over the whole 
bandwidth using a filter which has the inverse transfer function to that of the interstellar medium. 
The signal processing can be done on-line cither using finite impulse response filter devices 16 
or completely in software 290 1297) . Off-line data reduction, while disk-space limited, allows for 
more flexible analysis schemes |21| . 

The maximum time resolution obtainable via coherent dedispersion is the inverse of the total 
receiver bandwidth. The current state of the art is the detection |111) of features on nanosecond 
timescales in pulses from the 33-ms pulsar B0531+21 in the Crab nebula shown in Figure 1211 
Simple light travel-time arguments can be made to show that, in the absence of relativistic beaming 
effects |1(J4) . these incredibly bright bursts originate from regions less than 1 m in size. 
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Figure 21: A 120 (is window centred on a coherently-dedispersed giant pulse from the Crab pulsar 
showing high-intensity nanosecond bursts. Figure provided by Tim Hankins 1111}/ . 
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4.2 The timing model 

Ideally, in order to model the rotational behaviour of the neutron star, we require TOAs measured 
by an inertial observer. An observatory located on Earth experiences accelerations with respect to 
the neutron star due to the Earth's rotation and orbital motion around the Sun and is therefore not 
in an inertial frame. To a very good approximation, the solar system centre-of-mass (barycentre) 
can be regarded as an inertial frame. It is now standard practice |130) to transform the observed 
TOAs to this frame using a planetary ephemeris such as the JPL DE405 .293J- The transformation 
is summarised as the difference between barycentric (T) and observed (t) TOAs: 

T - t = ™ + ~ ^ + At Iei - At DM , (9) 
c led 

where r is the position of the observatory with respect to the barycentre, s is a unit vector in the 
direction of the pulsar at a distance d and c is the speed of light. The first term on the right hand 
side of Equation © is the light travel time from the observatory to the solar system barycentre. 
Incoming pulses from all but the nearest pulsars can be approximated by plane wavefronts. The 
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second term, which represents the delay due to spherical wavefronts, yields the parallax and hence 
d. This has so far only been measured for five nearby millisecond pulsars |276l 13201 11741 I287| . The 
term At le \ represents the Einstein and Shapiro corrections due to general relativistic time delays 
in the solar system |18) . Since measurements can be carried out at different observing frequencies 
with different dispersive delays, TOAs are generally referred to the equivalent time that would be 
observed at infinite frequency. This transformation is the term A^dm (see also Equation (|T|)). 

Following the accumulation of a number of TOAs, a surprisingly simple model is usually suffi- 
cient to account for the TOAs during the time span of the observations and to predict the arrival 
times of subsequent pulses. The model is a Taylor expansion of the rotational frequency Q — 2-k/P 
about a model value Slo at some reference epoch 7o. The model pulse phase is 

<f>(T) =<f> + (T- T )n + i(T - To) 2 f2 + . . • , (f 0) 

where T is the barycentric time and <^>o is the pulse phase at Tq. Based on this simple model, and 
using initial estimates of the position, dispersion measure and pulse period, a "timing residual" is 
calculated for each TOA as the difference between the observed and predicted pulse phases. 

A set of timing residuals for the nearby pulsar B1133+16 spanning almost 10 years is shown 
in Figure 1221 Ideally, the residuals should have a zero mean and be free from any systematic 
trends (see Panel a of Figure l2"2")l. To reach this point, however, the model needs to be refined in a 
bootstrap fashion. Early sets of residuals will exhibit a number of trends indicating a systematic 
error in one or more of the model parameters, or a parameter not incorporated into the model. 

From Equation (|10l) . an error in the assumed f2o results in a linear slope with time. A parabolic 
trend results from an error in (Iq (see Panel b of Figure Additional effects will arise if the 
assumed position of the pulsar (the unit vector § in Equation 0) used in the barycentric time 
calculation is incorrect. A position error results in an annual sinusoid (see Panel c of Figure |2"2"|) . 
A proper motion produces an annual sinusoid of linearly increasing magnitude (see Panel d of 
Figure EU 

After a number of iterations, and with the benefit of a modicum of experience, it is possible to 
identify and account for each of these various effects to produce a "timing solution" which is phase 
coherent over the whole data span. The resulting model parameters provide spin and astrometric in- 
formation with a precision which improves as the length of the data span increases. Timing observa- 
tions of the original millisecond pulsar B1937+21, spanning almost 9 years (exactly 165,711,423,279 
rotations!), measure a period of 1.5578064688197945 ± 0.0000000000000004 ms HE2IE21 defined 
at midnight UT on December 5 1988! Astrometric measurements are no less impressive, with 
position errors of ~ 20 /^arcsec being presently possible |333| . 

4.3 Timing stability 

Ideally, after correctly applying a timing model, we would expect a set of uncorrelated timing 
residuals with a zero mean and a Gaussian scatter with a standard deviation consistent with the 
measurement uncertainties. As can be seen in Figure l2"31 this is not always the case; the residuals 
of many pulsars exhibit a quasi-periodic wandering with time. 

Such "timing noise" is most prominent in the youngest of the normal pulsars |213l ITS] and 
present at a lower level in the much older millisecond pulsars |155lfTT) . While the physical processes 
of this phenomenon are not well understood, it seems likely that it may be connected to superfluid 
processes and temperature changes in the interior of the neutron star [3], or to processes in the 
magnetosphere jBHIEH]- 

The relative dearth of timing noise for the older pulsars is a very important finding. It implies 
that, presently, the measurement precision depends primarily on the particular hardware con- 
straints of the observing system. Consequently, a large effort in hardware development is presently 
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Figure 22: Timing model residuals for PSR B1133+16. Panel a: Residuals obtained from the 
best-fitting model which includes period, period derivative, position and proper motion. Panel b: 
Residuals obtained when the period derivative term is set to zero. Panel c: Residuals showing 
the effect of a 1-arcmin position error. Panel d: Residuals obtained neglecting the proper motion. 
The lines in Panels b-d show the expected behaviour in the timing residuals for each effect. Data 
provided by Andrew Lyne. 
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Figure 23: Examples of timing residuals for a number of normal pulsars. Note the varying scale 
on the ordinate axis, the pulsars being ranked in increasing order of timing "activity". Data taken 
from the Jodrell Bank timing program \2H1\ \124\ - Figure provided by Andrew Lyne. 
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being made to improve the precision of these observations using, in particular, coherent dedisper- 
sion outlined in Sect ion |4~T1 Much progress in this area has been made by groups at Princeton [258 , 
Berkeley Jodrell Bank UBC [323, Swinburne [H] and ATNF Q2]. From high qual- 

ity observations spanning over a decade [2741 12751 1155) . these groups have demonstrated that the 
timing stability of millisecond pulsars over such time-scales is comparable to terrestrial atomic 
clocks. 

This phenomenal stability is demonstrated in Figure l2*4l which shows a z 214 , a parameter 
closely resembling the Allan variance used by the clock community to estimate the stability of 
atomic clocks [3H51IT]. Both PSRs B1937+21 and B1855+09 seem to be limited by a power law 
component which produces a minimum in o~ z after 2 yr and 5 yr respectively. This is most likely 
a result of a small amount of intrinsic timing noise |155| . Although the baseline for the bright 
millisecond pulsar J0437— 4715 is shorter, its o~ z is already an order of magnitude smaller than the 
other two pulsars or the atomic clocks. Timing observations of an array of millisecond pulsars in the 
context of detecting gravitational waves from the Big Bang are discussed further in Section 14.5.31 
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Figure 24: The fractional stability of three millisecond pulsars compared to an atomic clock. Both 
PSRs B1855+09 and B1937+21 are comparable, or just slightly worse than, the atomic clock 
behaviour over timescales of a few years 12141- More recent timing of the millisecond pulsar 
304-31—4115 '1261 indicates that it is more stable than the atomic clock. 



4.4 Timing binary pulsars 

For binary pulsars, the simple timing model introduced in Section 14.21 needs to be extended to 
incorporate the additional motion of the pulsar as it orbits the common centre-of-mass of the 
binary system. Treating the binary orbit using Kepler's laws to refer the TOAs to the binary 
barycentre requires five additional model parameters: the orbital period P\>, projected semi-major 
orbital axis x, orbital eccentricity e, longitude of periastron to and the epoch of periastron passage 
T . This description, using five "Keplerian parameters", is identical to that used for spectroscopic 
binary stars. Analogous to the radial velocity curve in a spectroscopic binary, for binary pulsars 
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the orbit is described by the apparent pulse period against time. An example of this is shown in 
Panel a of Figure Alternatively, when radial accelerations can be measured, the orbit can also 
be visualised in a plot of acceleration versus period as shown in Panel b of Figure 12*51 This method 
is particularly useful for determining binary pulsar orbits from sparsely sampled data |102| . 
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Figure 25: Panel a: Keplerian orbital fit to the 669-day binary pulsar J04-07+1607 11901/ . Panel b: 
Orbital fit in the period- acceleration plane for the globular cluster pulsar 47 Tuc S \102L 

Constraints on the masses of the pulsar m p and the orbiting companion m c can be placed by 
combining x and Pb to obtain the mass function 
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where G is Newton's gravitational constant and i is the (initially unknown) angle between the 
orbital plane and the plane of the sky (i.e. an orbit viewed edge-on corresponds to i = 90°). In 
the absence of further information, the standard practice is to consider a random distribution of 
inclination angles. Since the probability that i is less than some value io is p(< zo) = 1 — cos(«o), 
the 90% confidence interval for i is 26° < i < 90°. For an assumed pulsar mass, the 90% confidence 
interval for m c can be obtained by solving Equation (fill) for i = 26° and 90°. If the sum of the 
masses M = m p + m c can be determined (e.g., through a measurement of relativistic periastron 
advance described below), then the condition sini < 1 sets a lower limit on the companion mass 
m c > (fmassM 2 ) 1 / 3 and a corresponding upper limit on the pulsar mass. 

Although most of the presently known binary pulsar systems can be adequately timed using 
Kepler's laws, there are a number which require an additional set of "post-Keplerian" (PK) pa- 
rameters which have a distinct functional form for a given relativistic theory of gravity. In general 
relativity (GR) the PK formalism gives the relativistic advance of periastron 
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the time dilation and gravitational redshift parameter 
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gj T^ /3 M- 4 / 3 m c (m p + 2m c ), (13) 
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and the two Shapiro delay parameters 

r = T Q m c 



(15) 



and 

^x^Y^Te^M^rn; 1 (16) 

which describe the delay in the pulses around superior conjunction where the pulsar radiation 
traverses the gravitational well of its companion. In the above expressions, all masses are in solar 
units, M = m p + m c , x = a p s'mi/c, s = smi and Tq = GM Q /c 3 — 4.925490947 fis. Some 
combinations, or all, of the PK parameters have now been measured for a number of binary pulsar 
systems. Further PK parameters due to aberration and relativistic deformation JSIJJ are not listed 
here but may soon be important for the double pulsar [185] . 

The key point in the PK definitions is that, given the precisely measured Kcplcrian parameters, 
the only two unknowns are the masses of the pulsar and its companion, m p and m c . Hence, from 
a measurement of just two PK parameters (e.g., lu and 7) one can solve for the two masses and, 
using Equation (|11|) . find the orbital inclination angle i. If three (or more) PK parameters are 
measured, the system is "overdetermined" and can be used to test GR (or, more generally, any 
other theory of gravity) by comparing the third PK parameter with the predicted value based on 
the masses determined from the other two. 

The first binary pulsar used to test GR in this way was PSR B1913+16 discovered by Hulsc 
& Taylor in 1974 129 . Measurements of three PK parameters (lu, 7 and J\>) were obtained from 
long-term timing observations at Arecibo |312l I313| . The measurement of orbital decay, which 
corresponds to a shrinkage of about 3.2 mm per orbit, is seen most dramatically as the gradually 
increasing shift in orbital phase for periastron passages with respect to a non-decaying orbit shown 
in Figure |2S| This figure includes recent Arecibo data taken since the upgrade of the telescope in 
the mid 1990s. The measurement of orbital decay, now spanning a 30-yr baseline [336 , is within 
0.2% of the GR prediction and provided the first indirect evidence for the existence of gravitational 
waves. Hulse and Taylor were awarded the 1993 Nobel Physics prize [3141 11281 1309) in recognition 
of their discovery of this remarkable laboratory for testing GR. 

More recently, five PK parameters have been measured for PSRs B1534+12 |298j and J0737— 
3039A |161) . For PSR B1534+12, the test of GR comes from measurements of lu, 7 and s, where 
the agreement between theory and observation is within 0.7% |298| . This test will improve in 
the future as the timing baseline extends and a more significant measurement of r can be made. 
Although a significant measurement of exists, it is known to be contaminated by kinematic 
effects which depend on the assumed distance to the pulsar 293 . Assuming GR to be correct, 
the observed and theoretical Pb values can be reconciled to provide a "relativistic measurement" 
of the distance d — 1.04 ± 0.03 kpc 288 . Prospects for independent parallax measurements of the 
distance to this pulsar using radio interferometry await more sensitive telescopes |289) . 

For PSR J0737— 3039, where two independent pulsar clocks can be timed, five PK parameters of 
the 22.7-ms pulsar "A" have been measured as well as two additional constraints from the measured 
mass function and projected semi-major axis of the 2.7-s pulsar "B". In terms of a laboratory for 
GR, then, J0737— 3039 promises to go well beyond the results possible from PSRs B1913+16 
and B1534+12. A useful means of summarising the limits so far is Figure 1771 which shows the 
allowed regions of parameter space in terms of the masses of the two pulsars. The shaded regions 
are excluded by the requirement that sinz < 1. Further constraints are shown as pairs of lines 
enclosing permitted regions as predicted by GR. The measurement of lu = 16.899 ±0.001 deg yr _1 
gives the total system mass M — 2.587 ± 0.001 Mq. The measurement of the projected semi-major 
axes of both orbits gives the mass ratio R = 1.071 ± 0.001. The mass ratio measurement is unique 
to the double pulsar system and rests on the basic assumption that momentum is conserved. This 
constraint should apply to any reasonable theory of gravity. The intersection between the lines for 
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Figure 26: Orbital decay in the binary pulsar B1913+16 system demonstrated as an increasing 
orbital phase shift for periastron passages with time. The GR prediction due entirely to the emission 
of gravitational radiation is shown by the parabola. Figure provided by Joel Weisberg '3361. 
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Figure 27: 'Mass-mass' diagram showing the observational constraints on the masses of the neutron 
stars in the double pulsar system J0737—3039. Inset is an enlarged view of the small square 
encompassing the intersection of the tightest constraints. Figure provided by Michael Kramer \lbl\f . 
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Co and R yield the masses of A and B as m A = 1.338± 0.001 M Q and m B = 1.249± 0.001 Mq. From 
these values, using Equations (fTSHTHf) the expected values of 7, Pb, r and s may be calculated and 
compared with the observed values. These four tests of GR all agree with the theory to within 
the uncertainties. Currently the tightest constraint is the Shapiro delay parameter s where the 
observed value is in agreement with GR at 0.1% level. 

Less than two years after its discovery, the double pulsar system has already surpassed the 
three decades of monitoring PSR B1913+16 and over a decade of timing PSR B1534+12 as a 
precision test of GR. On-going precision timing measurements of the double pulsar system should 
soon provide even more stringent and new tests of GR. Crucial to these measurements will be the 
timing of the 2.7-s pulsar B, where the observed profile is significantly affected by A's relativistic 
wind |198l 1216) . A careful decoupling of these profile variations is required to accurately measure 
TOAs for this pulsar and determine the extent to which the theory-independent mass ratio R can 
be measured. 

The relativistic effects observed in the double pulsar system are so large that corrections to 
higher post-Newtonian order may soon need to be considered. For example, Co may be measured 
precisely enough to require terms of second post-Newtonian order to be included in the compu- 
tations I n addition, in contrast to Newtonian physics, GR predicts that the spins of the 
neutron stars affect their orbital motion via spin-orbit coupling. This effect would most clearly 
be visible as a contribution to the observed Co in a secular |2fc>) and periodic fashion |337| . For 
the J0737— 3039 system, the expected contribution is about an order of magnitude larger than for 
PSR B1913+16 |198| . As the exact value depends on the pulsars' moment of inertia, a potential 
measurement of this effect allows the moment of inertia of a neutron star to be determined for the 
first time |81) . Such a measurement would be invaluable for studies of the neutron star equation 
of state and our understanding of matter at extreme pressure and densities |168) . 

The systems discussed above are all double neutron star binaries. A further self-consistency 
test of GR has recently been made in the 4.7-hr relativistic binary J1141— 6545, where the mea- 
surement 25^ of Co, 7 and Pb yield a pulsar mass of 1.30 ± 0.02 Mq and a companion mass 
of 0.99 ± 0.02 Mq. Since the mass of the companion is some seven standard deviations below 
the mean neutron star mass (see Figure I28f) . it is most likely a white dwarf. The observed 
i\ = ~(4±1) x 10~ 13 is consistent, albeit with limited precision, with the predicted value from GR 
(—3.8 x 10~ 13 ). Continued timing should reduce the relative error in Pb down to 1% by 2010 |25| . 

PK parameters have now been measured for a number of other binary pulsars which provide 
interesting constraints on neutron star masses |318l I292| . Figure [2H1 shows the distribution taken 
from a recent compilation 292 . While the young pulsars and the double neutron star binaries 
are consistent with, or just below, the canonical 1.4 Mq, we note that the millisecond pulsars in 
binary systems have, on average, significantly larger masses. This provides strong support for their 
formation through an extended period of accretion in the past, as discussed in Section l2~6l 

4.5 Pulsar timing and gravitational wave detection 

We have seen in the previous Sectior fOl how pulsar timing can be used to provide indirect evidence 
for the existence of gravitational waves from coalescing stellar-mass binaries. In this final section, 
we look at how pulsar timing might soon be used to detect gravitational radiation directly. The idea 
to use pulsars as natural gravitational wave detectors was first explored in the late 1970s [2771 155] , 
The basic concept is to treat the solar system barycentre and a distant pulsar as opposite ends of 
an imaginary arm in space. The pulsar acts as the reference clock at one end of the arm sending out 
regular signals which are monitored by an observer on the Earth over some time-scale T. The effect 
of a passing gravitational wave would be to perturb the local spacetime metric and cause a change 
in the observed rotational frequency of the pulsar. For regular pulsar timing observations with 
typical TO A uncertainties of ctoa, this "detector" would be sensitive to waves with dimcnsionless 
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Figure 28: Distribution of neutron star masses as inferred from timing observations of binary 
pulsars '2921. The vertical dotted line shows the canonical neutron star mass of IAMq. Figure 
provided by Ingrid Stairs. 
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amplitudes h > etoa / T and frequencies / ~ 1/T I39| . 
4.5.1 Probing the gravitational wave background 

Many cosmological models predict that the Universe is presently filled with a low-frequency stochas- 
tic gravitational wave background (GWB) produced during the big bang era |243| . A significant 
component |2fiOI I134| is the gravitational radiation from massive black hole mergers. In the ideal 
case, the change in the observed frequency caused by the GWB should be detectable in the set of 
timing residuals after the application of an appropriate model for the rotational, astrometric and, 
where necessary, binary parameters of the pulsar. As discussed in Section l4.2l all other effects be- 
ing negligible, the rms scatter of these residuals a would be due to the measurement uncertainties 
and intrinsic timing noise from the neutron star. 

For a GWB with a flat energy spectrum in the frequency band / ± //2 there is an additional 
contribution to the timing residuals cr g |85) . When fT 3> 1, the corresponding wave energy density 
is 

243tt 3 / 4 ^ 

An upper limit to p g can be obtained from a set of timing residuals by assuming the rms scatter 
is entirely due to this effect (er = er g ). These limits are commonly expressed as a fraction of the 
energy density required to close the Universe 

Pc = H ~ 2 x 1(T 29 /, 2 g cm" 3 , (18) 

where the Hubble constant Hq = 100 h km s _1 Mpc. 

This technique was first applied |273) to a set of TOAs for PSR B1237+25 obtained from 
regular observations over a period of 11 years as part of the JPL pulsar timing programme |89j . 
This pulsar was chosen on the basis of its relatively low level of timing activity by comparison with 
the youngest pulsars, whose residuals are ultimately plagued by timing noise (see Section I4.3J) . By 
ascribing the rms scatter in the residuals (a — 240 ms) to the GWB, the limit p g /p c ^ 4x 10~ 3 h~ 2 
for a centre frequency / = 7 nHz. 

This limit, already well below the energy density required to close the Universe, was further 
reduced following the long-term timing measurements of millisecond pulsars at Arecibo (see Sec- 
tion I n the intervening period, more elaborate techniques had been devised (311 1391 T303I to 
look for the likely signature of a GWB in the frequency spectrum of the timing residuals and to 
address the possibility of "fitting out" the signal in the TOAs. Following |SJ it is convenient to 
define 

1 dlogpg 
Pc dlog/ 

i.e. the energy density of the GWB per logarithmic frequency interval relative to p c . With this 
definition, the power spectrum of the GWB [1251 13*9"] is 



" 6 = -S (19) 



where / yr -i is frequency in cycles per year. The long-term timing stability of B1937+21, discussed 
in Section l4~3l limits its use for periods > 2 yr. Using the more stable residuals for PSR B1855+09, 
and a rigorous statistical analysis |319) . the current 95% confidence upper limit is Vl g h 2 < 10 -8 . 
This limit is difficult to reconcile with most cosmic string models for galaxy formation |45l I319j . 

For binary pulsars, the orbital period provides an additional clock for measuring the effects of 
gravitational waves. In this case, the range of frequencies is not limited by the time span of the 
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observations, allowing the detection of waves with periods as large as the light travel time to the 
binary system j^]. The most stringent results presently available are based on the B1855+09 limit 
n g h 2 < 2.7 x 1CT 4 in the frequency range KT 11 < / < 4.4 x 10" 9 Hz [ToD] . 

4.5.2 Constraints on massive black hole binaries 

In addition to probing the GWB, pulsar timing is beginning to place interesting constraints on 
the existence of massive black hole binaries. Arecibo data for PSRs B1937+21 and J1713+0747 
already make the existence of an equal-mass black hole binary in Sagittarius A* unlikely |176) . 
More recently, timing data from B1855+09 have been used to virtually rule out the existence of a 
proposed supermassive black hole as the explanation for the periodic motion seen at the centre of 
the radio galaxy 3C66B pM] . 

A simulation of the expected modulations of the timing residuals for the putative binary system, 
with a total mass of 5.4 x 10 10 Mq, is shown along with the observed timing residuals in Figure 1291 
Although the exact signature depends on the orientation and eccentricity of the binary system, 
Monte Carlo simulations show that the existence of such a massive black hole binary is ruled out 
with at least 95% confidence |135j . 
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Figure 29: Top panel: Observed timing residuals for PSR B1855+09. Bottom panel: Simulated 
timing residuals induced from a putative black hole binary in 3C66B. Figure provided by Rick 
Jenet f7%f . 



4.5.3 A millisecond pulsar timing array 

A natural extension of the single-arm detector concept discussed above is the idea of using timing 
data for a number of pulsars distributed over the whole sky to detect gravitational waves |118| . Such 
a "timing array" would have the advantage over a single arm in that, through a cross-correlation 
analysis of the residuals for pairs of pulsars distributed over the sky, it should be possible to 
separate the timing noise of each pulsar from the signature of the GWB common to all pulsars in 
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the array. To illustrate this, consider the fractional frequency shift of the ith pulsar in an array 



Vi 



(21) 



where on is a geometric factor dependent on the line-of-sight direction to the pulsar and the 
propagation and polarisation vectors of the gravitational wave of dimensionless amplitude A. The 
timing noise intrinsic to the pulsar is characterised by the function Mi. The result of a cross- 
correlation between pulsars i and j is then 



ij(A 2 ) + ai(AAfj) + a J (AM l ) + (MiMj), 



(22) 



where the bracketed terms indicate cross-correlations. Since the wave function and the noise 
contributions from the two pulsars are independent, the cross correlation tends to ctiaj(A 2 ) as the 
number of residuals becomes large. Summing the cross-correlation functions over a large number of 
pulsar pairs provides additional information on this term as a function of the angle on the sky |117| 
and allows the separation of the effects of clock and ephemeris errors from the GWB |97) . 
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Figure 30: Summary of the gravitational wave spectrum showing the location in phase space of 
the pulsar timing array (PTA) and its extension with the Square Kilometre Array (SKA). Figure 
updated by Michael Kramer ]16S\/ from an original design by Richard Battye. 
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A preliminary analysis applying the timing array concept to Arecibo data for three millisecond 
pulsars (B1937+21, B1855+09 and J1737+0747) now spanning 17 years has reduced the energy 
density limit to fl g h 2 < 2 x 1CP 9 |175) 3 . This limit, and the region of the gravitational wave energy 
density spectrum probed by the current pulsar timing array is shown in Figure |3"U1 where it can be 
seen that the pulsar timing regime is complementary to the higher frequency bands of LISA and 
LIGO. 

A number of long-term timing projects are now underway to make a large-scale pulsar timing 
array a reality. At Arecibo and Green Bank, regular timing of a dozen or more millisecond pulsars 
are now carried out using on-line data acquisition systems by groups at Berkeley 325.J, Prince- 
ton |258) and UBC |323| . A similar system is also being commissioned at Jodrell Bank |142) . The 
Berkeley group have installed identical sets of data taking equipment at Effelsberg and Nangay 
which has enabled these telescopes to make regular high-quality millisecond pulsar timing obser- 
vations [1741 11)5] . In the southern hemisphere, a high-precision timing program at Parkes 333 
has been in existence since the discovery bright nearby millisecond pusar J0437— 4715 |138) . Since 
February 2004, weekly Parkes observations of ^20 millisecond pulsars are already achieving a limit 
of VL g h 2 < 10~ 4 , with the ultimate goal being to reach a limits of VL g h 2 < 10~ 10 before 2010 |121) . 

Looking further ahead, the increase in sensitivity provided by the Square Kilometre Array [1321 
1162) should further improve the limits of the spectrum probed by pulsar timing. As Figure l3"Ul shows . 
the SKA could provide up to two orders of magnitude improvement over current capabilities. 

4.6 Going further 

For further details on the techniques and prospects of pulsar timing, a number of excellent reviews 
are available |18l 13081 1281 12911 1185) . Two audio files and slides from lectures at the centennial 
meeting of the American Physical Society are also available on-line [151 1341) . The tests of general 
relativity discussed in Section [4.41 along with further effects, are reviewed elsewhere in this journal 
by Will [342) and Stairs |291) . Further discussions on the realities of using pulsars as gravity wave 
detectors can be found in other review articles |271l[THIl21| . 



3 It is likely, however, that a more detailed statistical analysis of the data are required to confirm or refute this 
result E2I- 
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5 Summary and Future Prospects 



The main aim of this article was to review some of the many astrophysical applications provided 
by the present sample of binary and millisecond radio pulsars. The topics covered here, along 
with the bibliography and associated tables of observational parameters, should be useful to those 
wishing to delve deeper into the vast body of literature that exists. 

Through an understanding of the Galactic population of radio pulsars summarised in Section |21 
it is possible to predict the detection statistics of terrestrial gravitational wave detectors to nearby 
rapidly spinning neutron stars (see Section [3.3[l . as well as coalescing relativistic binaries at cosmic 
distances (see Section l3.4[) . Continued improvements in gravitational wave detector sensitivities 
should result in a number of interesting developments and contributions in this area. These devel- 
opments and contributions might include the detection of presently known radio pulsars, as well 
as a population of coalescing binary systems which have not yet been detected as radio pulsars. 

The phenomenal timing stability of radio pulsars leads naturally to a large number of appli- 
cations, including their use as laboratories for relativistic gravity (see Section WM and as natural 
detectors of gravitational radiation (see Section [4.5(1 . Long-term timing experiments of the present 
sample of millisecond and binary pulsars currently underway appear to have tremendous potential 
in these areas and perhaps detect the gravitational wave background (if it exists) within the next 
decade. 

These applications will benefit greatly from the continuing discovery of new systems by the 
present generation of radio pulsar searches which continue to probe new areas of parameter space. 
Based on the results presented in Section l3~3l it is clear that we are aware of only about 1% of 
the total active pulsar population in our Galaxy. It is therefore likely that we have not seen all of 
the pulsar zoo. More sensitive surveys using a multibeam system on the Arecibo telescope |227] 
are now being carried out. Future surveys with the Square Kilometre Array |132| will ultimately 
provide a far more complete census of the Galactic pulsar population. With the double pulsar 
now checked off the list of pulsar milestones given in the last version of this review 180 , we look 
forward to new discoveries. Two possible "holy grails" of pulsar astronomy which may soon be 
found are the following ones: 

A pulsar— black hole binary system: Following the wide variety of science possible from the 
new double pulsar J0737— 3039 (see Sections 12.61 13.4. Tl and l4~4")) . a radio pulsar with a black- 
hole companion would without doubt be a fantastic gravitational laboratory. Excellent places 
to look for such systems are globular clusters and the Galactic centre [247] . 

A sub-millisecond pulsar: With the new discoveries in Terzan 5, it now appears that the 
original millisecond pulsar, B1937+21 is no longer the most rapidly rotating neutron star 
known |233j . Do neutron stars with kHz spin rates exist? Searches now have sensitivity to 
such objects |41] and a discovery of even one would constrain the equation of state of matter 
at high densities. As mentioned in Section 12.21 the R-mode instability may prevent neutron 
stars from reaching such high spin rates [361 1601 17). 

Pulsar astronomy remains an extremely active area of modern astrophysics and the next decade 
will undoubtedly continue to produce new results from currently known objects as well as new 
surprises. In my opinion, pulsar research is currently limited by a shortage of researchers, and not 
necessarily resources. Keen graduate students more than ever are needed to help shape the future 
of this exciting and continually evolving field. 
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A Tables of Binary and Millisecond Pulsars 
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Tabic 2: Parameters for the 15 isolated millisecond pulsars currently known in the Galactic disk. 
Two of the pulsars (30609+2130 and J2235+1506) are thought to be the descendants of high-mass 
binary systems. Listed are the spin period P, the base-10 logarithms of the characteristic age r c 
and surface magnetic field strength B (see Section \2.2ty . the distance d derived from the Cordes & 
Lasio electron density model \73\ \ 74\j or independently (when available), and the transverse speed 
v t inferred from d and a proper motion measurement (when available). Key publications for each 
pulsar are referenced to the bibliography. 
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Table 3: Parameters for the 15 eccentric (e > 0.05J binary pulsars currently known that are not in globular clusters. In addition to the 
parameters listed in Tabled we also give the binary period P^, the projected semi-major axis of the orbit x in units of light seconds, the 
orbital eccentricity e, and the companion mass mi evaluated from the mass function assuming a pulsar mass o/1.4Mq and an inclination 
angle of 60 degrees (see Section \4-4{ ) or (when known) from independent measurements. 
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Table 4: Parameters for 50 low- eccentricity binary pulsars currently known in the Galactic disk. 
The parameters listed are defined in Tables\^ and\^ 
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Table 5: Parameters for the 108 pulsars currently known in globular clusters. In addition to the 
parameters defined in Tables\^ and\^ we also list the measured spin period derivative P and name 
of associated cluster. 



References 



[1] Allan, D.W., "The Allan Variance", personal homepage, Allan's Time Interval Metrology 
Enterprises, (1999). URL (cited on 7 December 2000): 
: //www . allanstime . com/AllanVar iance/ 

[2] Alpar, M.A., Cheng, A.F., Ruderman, M.A., and Shaham, J., "A new class of radio pulsars", 
Nature, 300, 728-730, (1982). 

[3] Alpar, M.A., Nandkumar, R., and Pines, D., "Vortex Creep and the Internal Temperature 
of Neutron Stars: Timing Noise of Pulsars", Astrophys. J., 311, 197-213, (1986). 

[4] Anderson, S.B., A study of recycled pulsars in globular clusters, Ph.D. Thesis, (California 
Institute of Technology, Pasadena, U.S.A., 1992). 

[5] Anderson, S.B., Gorham, P.W., Kulkarni, S.R., Prince, T.A., and Wolszczan, A., "Discovery 
of two radio pulsars in the globular cluster M15", Nature, 346, 42-44, (1990). 

[6] Anderson, S.B., Wolszczan, A., Kulkarni, S.R., and Prince, T.A., "Observations of Two 
Millisecond Pulsars in the Globular Cluster NGC 5904" , Astrophys. J., 482, 870-873, (1997). 

[7] Andersson, N., Jones, D.I., Kokkotas, K., and Stergioulas, N., "r-Mode runaway and rapidly 
rotating neutron stars", Astrophys. J. Lett., 534, L75-L78, (2000). 

[8] Arzoumanian, Z., Chcrnoff, D. F., and Cordes, J. M., "The Velocity Distribution of Isolated 
Radio Pulsars", Astrophys. J., 568, 289-301, (2002). 

[9] Arzoumanian, Z., Cordes, J.M., and Wasserman, I., "Pulsar Spin Evolution, Kinematics, 
and the Birthrate of Neutron Star Binaries", Astrophys. J., 520, 696-705, (1999). 

[10] Arzoumanian, Z., Fruchter, A.S., and Taylor, J.H., "Orbital Variability in the Eclipsing 
Pulsar Binary PSR B1957+20", Astrophys. J. Lett, 426, L85-L88, (1994). 

[11] Arzoumanian, Z., Nice, D.J., Taylor, J.H., and Thorsett, S.E., "Timing Behavior of 96 Radio 
Pulsars", Astrophys. J., 422, 671-680, (1994). 

[12] Australia Telescope National Facility, "ATNF Pulsar Group", project homepage, (2005). 

URL (cited on 22 January 2005): 

http : //www . atnf . csiro . au/research/pulsar 

[13] Australia Telescope National Facility, "The Pulsar Catalogue" , web interface to database, 

(2005). URL (cited on 17 January 2005): 

http : //www . atnf . csiro . au/research/pulsar/psrcat 

[14] Backer, D.C., "Timing of Millisecond Pulsars", in van Paradijs, J., van del Heuvel, E.P.J., 
and Kuulkers, E., eds., Compact Stars in Binaries (IAU Symposium 165), Proceedings of the 
165th Symposium of the International Astronomical Union, held in the Hague, The Nether- 
lands, August 15-19, 1994, 197-211, (Kluwer, Dordrecht, Netherlands; Boston, U.S.A., 
1996). 

[15] Backer, D.C., "Pulsars - Nature's Best Clocks", lecture notes, American Physical Society, 

(1999). URL (cited on 8 December 2000): 

http : //www . apscenttalks . org/pres_masterpage . cf m?nameID=65 



54 



[16] Backer, D.C., Dexter, M.R., Zepka, A.F., Ng, D., Werthimcr, D.J., Ray, P.S., and Foster, 
R.S., "A Programmable 36 MHz Digital Filter Bank for Radio Science", Publ. Astron. Soc 
Pac, 109, 61-68, (1997). 

[17] Backer, D.C., Foster, R.F., and Sallmcn, S., "A second companion of the millisecond pulsar 
1620-26", Nature, 365, 817-819, (1993). 

[18] Backer, D.C., and Hellings, R.W., "Pulsar Timing and General Relativity", Annu. Rev. 
Astron. Astrophys., 24, 537-575, (1986). 

[19] Backer, D.C., Kulkarni, S.R., Hciles, C, Davis, M.M., and Goss, W.M., "A millisecond 
pulsar", Nature, 300, 615-618, (1982). 

[20] Bailes, M., "The Origin of pulsar velocities and the velocity-magnetic moment correlation.", 
Astrophys. J., 342, 917-927, (1989). 

[21] Bailes, M., "Precision timing at the Parkes 64-m radio telescope", in Bailes, M., Nice, D.J., 
and Thorsett, S.E., eds., Radio Pulsars, Proceedings of meeting held at Mediterranean Ago- 
nomic Institute of Chania, Crete, Greece, 26-29 August 2002, vol. 302 of ASP Conference 
Series, 57-64, (Astronomical Society of the Pacific, San Fransisco, U.S.A., 2003). 

[22] Bailes, M., Harrison, P.A., Lorimer, D.R., Johnston, S., Lync, A.G., Manchester, R.N., 
DAmico, N., Nicastro, L., Tauris, T.M., and Robinson, C, "Discovery of Three Galactic 
Binary Millisecond Pulsars", Astrophys. J. Lett., 425, L41-L44, (1994). 

[23] Bailes, M., Johnston, S., Bell, J.F., Lorimer, D.R., Stappers, B.W., Manchester, R.N., Lyne, 
A.G., D'Amico, N., and Gaensler, B.M., "Discovery of four isolated millisecond pulsars", 
Astrophys. J., 481, 386-391, (1997). 

[24] Bailes, M., Manchester, R.N., Kesteven, M.J., Norris, R.P., and Reynolds, J.E., "The proper 
motion of six southern radio pulsars", Mon. Not. R. Astron. Soc, 247, 322-326, (1990). 

[25] Bailes, M., Ord, S.M., Knight, H.S., and Hotan, A.W., "Self-Consistency of Rclativistic 
Observables with General Relativity in the White Dwarf-Neutron Star Binary PSR J1141- 
6545", Astrophys. J. Lett, 595, L49-L52, (2003). 

[26] Barker, B.M., and O'Connell, R.F., "Gravitational two-body problem with arbitrary masses, 
spins, and quadrupole moments", Phys. Rev. D, 12, 329-335, (1975). 

[27] Bclczynski, K., Kalogera, V.M., and Bulik, T., "A Comprehensive Study of Binary Com- 
pact Objects as Gravitational Wave Sources: Evolutionary Channels, Rates, and Physical 
Properties", Astrophys. J., 572, 407-431, (2002). 

[28] Bell, J.F., "Radio Pulsar Timing", Adv. Space Res., 21, 131-147, (1998). 

[29] Bell, J.F., Bailes, M., Manchester, R.N., Lyne, A.G., Camilo, F., and Sandhu, J.S., "Timing 
measurements and their implications for four binary millisecond pulsars", Mon. Not. R. 
Astron. Soc., 286, 463-469, (1997). 

[30] Bell, J.F., Bessell, M.S., Stappers, B.W., Bailes, M., and Kaspi, V.M., "PSR J0045-7319: A 
dual-line binary radio pulsar", Astrophys. J. Lett., 447, L117-L119, (1995). 

[31] Bcrtotti, B., Carr, B.J., and Rees, M.J., "Limits from the timing of pulsars on the cosmic 
gravitational wave background", Mon. Not. R. Astron. Soc, 203, 945-954, (1983). 



55 



[32] Beskin, V.S., Gurevich, A.V., and Istomin, Y.N., Physics of the Pulsar Magnetosphere, 
(Cambridge University Press, Cambridge, U.K.; New York, U.S.A., 1993). 

[33] Bhattacharya, D., and van den Heuvel, E.P.J. , "Formation and evolution of binary and 
millisecond radio pulsars", Phys. Rep., 203, 1-124, (1991). 

[34] Biggs, J.D., "Meridional compression of radio pulsar beams", Mon. Not. R. Astron. Soc., 
245, 514-521, (1990). 

[35] Biggs, J.D., Bailes, M., Lyne, A.G., Goss, W.M., and Fruchter, A.S., "Two radio pulsars in 
the globular cluster NGC 6624", Mon. Not. R. Astron. Soc., 267, 125-128, (1994). 

[36] Bildsten, L., "Gravitational radiation and rotation of accreting neutron stars", Astrophys. J. 
Lett, 501, L89-L93, (1998). 

[37] Bisnovatyi-Kogan, G.S., and Komberg, B.V., "Pulsars and Close Binary Systems", Sov. 
Astron., 18, 217-221, (1974). 

[38] Blaauw, A., "On the origin of the O- and B-type stars with high velocities (the 'Run-away' 
stars), and related problems", Bull. Astron. Inst. Neth., 15, 265-290, (1961). 

[39] Blandford, R.D., Narayan, R., and Romani, R.W., "Arrival-time analysis for a millisecond 
pulsar", J. Astrophys. Astron., 5, 369-388, (1984). 

[40] Boriakoff, V., Buccheri, R., and Fauci, F., "Discovery of a 6.1-ms binary pulsar, PSR 
1953+29", Nature, 304, 417-419, (1983). 

[41] Burderi, L., and D'Amico, N., "Probing the Equation of State of Ultradense Matter with a 
Submillisecond Pulsar Search Experiment", Astrophys. J., 490, 343-352, (1997). 

[42] Burgay, M., The Parkes High-Latitude Pulsar Survey and the Discovery of the first Double 
Pulsar, Ph.D. Thesis, (Bologna University, Bologna, Italy, 2004). 

[43] Burgay, M., Burderi, L., Possenti, A., D'Amico, N., Manchester, R.N., Lyne, A.G., Camilo, 
F., and Campana, S., "A Search for Pulsars in Quiescent Soft X-Ray Transients. I.", Astro- 
phys. J., 589, 902-910, (2003). 

[44] Burgay, M., D'Amico, N., Possenti, A., Manchester, R.N., Lyne, A.G., Joshi, B.C., McLaugh- 
lin, M.A., Kramer, M., Sarkissian, J.M., Camilo, F., Kalogera, V., Kim, C, and Lorimer, 
D.R., "An increased estimate of the merger rate of double neutron stars from observations 
of a highly relativistic system", Nature, 426, 531-533, (2003). 

[45] Caldwell, R.R., and Allen, B., "Cosmological constraints on cosmic-string gravitational ra- 
diation", Phys. Rev. D, 45, 3447-3468, (1992). 

[46] Caltech/MIT, "Laser Interferometric Gravitational Wave Observatory", project homepage, 
(2005). URL (cited on 29 January 2005): 
http : //www . ligo . caltech . edu 

[47] Camilo, F., "Millisecond Pulsar Searches", in Alpar, A., Kiziloglu, U., and van Paradis J., 
eds., The Lives of the Neutron Stars, Proceedings of the NATO Advanced Study Institute, 
Kemer, Turkey, August 29 - September 12, 1993, vol. 450 of NATO ASI Series, Series C, 
243-257, (Kluwer, Dordrecht, Netherlands; Boston, U.S.A., 1995). 

[48] Camilo, F., A Search for Millisecond Pulsars, Ph.D. Thesis, (Princeton University, Princeton, 
U.S.A., 1995). 



56 



[49] Camilo, F., "Intermediate-Mass Binary Pulsars: a New Class of Objects?", in Johnston, S., 
Walker, M.A., and Bailes, M., eds., Pulsars: Problems and Progress (IAU Colloquium 160), 
Proceedings of the 160th Colloquium of the IAU, held at the Research Centre for Theoretical 
Astrophysics, University of Sydney, Australia, 8-12 January 1996, vol. 105 of ASP Conference 
Series, 539-540, (Astronomical Society of the Pacific, San Francisco, U.S.A., 1996). 

[50] Camilo, F., "Present and future pulsar searches", in Jackson, N., and Davis, R.J., eds., 
High- Sensitivity Radio Astronomy, Proceedings of a meeting held at Jodrcll Bank, Univer- 
sity of Manchester, January 22-26, 1996, Cambridge Contemporary Astrophysics, 14-22, 
(Cambridge University Press, Cambridge, U.K.; New York, U.S.A., 1997). 

[51] Camilo, F., "Pulsar Searches in the Northern Hemisphere", in Arzoumanian, Z., van der 
Hooft, F., and van den Heuvel, E.P.J. , eds., Pulsar Timing, General Relativity, and the In- 
ternal Structure of Neutron Stars, Proceedings of the colloquium, Amsterdam, 24-27 Septem- 
ber 1996, 115, (Royal Netherlands Academy of Arts and Sciences, Amsterdam, Netherlands, 
1999). 

[52] Camilo, F., Foster, R.S., and Wolszczan, A., "High-precision timing of PSR J1713+0747: 
Shapiro delay", Astrophys. J. Lett, 437, L39-L42, (1994). 

[53] Camilo, F., Lorimer, D.R., Freire, P.C.C., Lyne, A.G., and Manchester, R.N., "Observations 
of 20 Millisecond Pulsars in 47 Tucanae at 20 Centimeters", Astrophys. J., 535, 975-990, 
(2000). 

[54] Camilo, F., Lyne, A.G., Manchester, R.N., Bell, J.F., Stairs, I.H., D'Amico, N., Kaspi, V.M., 
Posscnti, A., Crawford, F., and McKay, N.P.F., "Discovery of Five Binary Radio Pulsars", 
Astrophys. J. Lett., 548, L187-L191, (2001). 

[55] Camilo, F., Nice, D.J., and Taylor, J.H., "Discovery of two fast-rotating pulsars", Astrophys. 
J. Lett, 412, L37-L40, (1993). 

[56] Camilo, F., Nice, D.J., and Taylor, J.H., "A search for millisecond pulsars at galactic latitudes 
-50° < b < -20°", Astrophys. J., 461, 812-819, (1996). 

[57] Camilo, F., Nice, D.J., Taylor, J.H., and Shrauner, J. A., "Princeton-Arecibo Declination- 
Strip Survey for Millisecond Pulsars. L", Astrophys. J., 469, 819-827, (1996). 

[58] Camilo, F., and Rasio, F.A., "Radio Pulsars in Globular Clusters" , in Rasio, F.A., and Stairs, 
I.H., eds., Binary Radio Pulsars, Proceedings of meeting held at the Aspen Center for Physics, 
USA, 12 Janaury - 16 January 2004, vol. 328 of ASP Conference Series, (Astronomical Society 
of the Pacific, San Francisco, U.S.A., 2005). 

[59] Chakrabarty, D., and Morgan, E.H., "The 2 hour orbit of a binary millisecond X-ray pulsar" , 
Nature, 394, 346-348, (1998). 

[60] Chakrabarty, D., Morgan, E.H., Muno, M.P., Galloway, D.K., Wijnands, R., van der Klis, 
M., and Markwardt, C.B., "Nuclear-powered millisecond pulsars and the maximum spin 
frequency of neutron stars", Nature, 424, 42-44, (2003). 

[61] Champion, D.J., "Latest timing results for PSR J1829+2456", personal communication. 

[62] Champion, D.J., Lorimer, D.R., McLaughlin, M.A., Arzoumanian, Z., Cordes, J.M., Weis- 
berg, J., and Taylor, J.H., "PSR J1829+2456: a new relativistic binary system", Mon. Not. 
R. Astron. Soc., 350, L61-L65, (2004). 



57 



[63] Chandler, A.M., Pulsar Searches: From Radio to Gamma-Rays, Ph.D. Thesis, (California 
Institute of Technology, Pasadena, U.S.A., 2003). 

[64] Chen, K., and Ruderman, M., "Pulsar death rates and death valley", Astrophys. J., 402, 
264-270, (1993). 

[65] Cheng, K.S., "Could Glitches Inducing Magnetospheric Fluctuations Produce Low Frequency 
Timing Noise?", Astrophys. J., 321, 805-812, (1987). 

[66] Cheng, K.S., "Outer Magnetospheric Fluctuations and Pulsar Timing Noise", Astrophys. J., 
321, 799-804, (1987). 

[67] Clifton, T.R., Lyne, A.G., Jones, A.W., McKenna, J., and Ashworth, M., "A high frequency 
survey of the Galactic plane for young and distant pulsars", Mon. Not. R. Astron. Soc., 254, 
177-184, (1992). 

[68] Cognard, I., and Backer, D.C., "A Microglitch in the Millisecond Pulsar PSR B1821-24 in 
M28", Astrophys. J. Lett, 612, L125-L127, (2004). 

[69] Cooray, A., "Gravitational wave background of neutron star-white dwarf binaries" , Mon. 
Not. R. Astron. Soc., 354, 25-30, (2005). 

[70] Cordes, J.M., and Chernoff, D.F., "Neutron Star Population Dynamics. I. Millisecond Pul- 
sars", Astrophys. J., 482, 971-992, (1997). 

[71] Cordes, J.M., and Chernoff, D.F., "Neutron Star Population Dynamics. II. Three-dimensional 
Space Velocities of Young Pulsars", Astrophys. J., 505, 315-338, (1998). 

[72] Cordes, J.M., and Helfand, D.J., "Pulsar timing III. Timing Noise of 50 Pulsars" , Astrophys. 
J. Lett, 239, 640-650, (1980). 

[73] Cordes, J.M., and Lazio, T.J.W., "NE2001. I. A New Model for the Galactic Distribution of 
Free Electrons and its Fluctuations", (2002). URL (cited on 22 January 2005): 
http : //arXiv . org/abs/astro-ph/0207156 

[74] Cordes, J.M., and Lazio, T.J.W., "NE2001. II. Using Radio Propagation Data to Construct 
a Model for the Galactic Distribution of Free Electrons", (2003). URL (cited on 22 January 

2005): 

http : //arXiv . org/abs/astro-ph/0301598 

[75] Cornell University, "arXiv.org e-Print archive" , web interface to database, (2001) . URL (cited 
on 19 March 2001): 
http : //arXiv .orgl 

[76] D'Amico, N., Bailes, M., Lyne, A.G., Manchester, R.N., Johnston, S., Fruchter, A.S., and 
Goss, W.M., "PSR B1802-07: A globular cluster pulsar in an eccentric binary system", Mon. 
Not. R. Astron. Soc, 260, L7-L10, (1993). 

[77] D'Amico, N., Lyne, A.G., Manchester, R.N., Possenti, A., and Camilo, F., "Discovery of 
Short-Period Binary Millisecond Pulsars in Four Globular Clusters", Astrophys. J. Lett., 
548, L171-L174, (2001). 

[78] D'Amico, N., Possenti, A., Fici, L., Manchester, R.N., Lyne, A.G., Camilo, F., and 
Sarkissian, J., "Timing of Millisecond Pulsars in NGC 6752: Evidence for a High Mass- 
to-Light Ratio in the Cluster Core", Astrophys. J. Lett., 570, L89-L92, (2002). 



58 



[79] D'Amico, N., Possenti, A., Manchester, R.N., Sarkissian, J., Lyne, A.G., and Camilo, F., 
"An Eclipsing Millisecond Pulsar with a Possible Main-Sequence Companion in NGC 6397" , 
Astrophys. J. Lett, 561, L89-L92, (2001). 

[80] Damour, T., and Deruellc, N., "General Relativistic Celestial Mechanics of Binary Systems. 
II. The Post-Newtonian Timing Formula" , Ann. Inst. Henri Poincare A, 44, 263-292, (1986). 

[81] Damour, T., and Schafer, G., "Higher-Order Relativistic Periastron Advances and Binary 
Pulsars", Nuovo Cimento, 101, 127, (1988). 

[82] Damour, T., and Vilenkin, A., "Gravitational radiation from cosmic (super)strings: Bursts, 
stochastic background, and observational windows" , Phys. Rev. D, 71, 063510-1-13, (2005). 
Related online version (cited on 30 January 2005): 
http : //arXiv . org/abs/hep-th/0410222 

[83] Dana, P.H., "The Global Positioning System", project homepage, Department of Geography, 
University of Colorado at Boulder, (2000). URL (cited on 19 March 2001): 
http : //www . Colorado . edu/geography/gcraf t/notes/gps/gps_f . html 

[84] Davis, M.M., Taylor, J.H., Weisberg, J.M., and Backer, D.C., "High-precision timing obser- 
vations of the millisecond pulsar PSR 1937+21", Nature, 315, 547-550, (1985). 

[85] Detweiler, S.L., "Pulsar timing measurements and the search for gravitational waves", As- 
trophys. J., 234, 1100-1104, (1979). 

[86] Dewey, R.J., and Cordes, J.M., "Monte Carlo Simulations of Radio Pulsars and their Pro- 
genitors", Astrophys. J., 321, 780-798, (1987). 

[87] Dewey, R.J., Maguire, CM., Rawley, L.A., Stokes, G.H., and Taylor, J.H., "Binary pulsar 
with a very small mass function", Nature, 322, 712-714, (1986). 

[88] Dewey, R.J., Stokes, G.H., Segelstein, D.J., Taylor, J.H., and Weisberg, J.M., "The Period 
Distribution of Pulsars", in Reynolds, S.P., and Stinebring, D.R., eds., Birth and Evolution 
of Neutron Stars: Issues Raised by Millisecond Pulsars, Proceedings of a workshop, held at 
the National Radio Astronomy Observatory, Green Bank, West Virginia, June 6-8, 1984, 
vol. 8 of NRAO Workshop, 234-240, (National Radio Astronomy Observatory, Green Bank, 
U.S.A., 1984). 

[89] Downs, G.S., and Reichley, P.E., "JPL pulsar timing observations. II. Geocentric arrival 
times", Astrophys. J. Suppl. Ser., 53, 169-240, (1983). 

[90] Edwards, R., and Bailes, M., "Discovery of two relativistic neutron star- white dwarf bina- 
ries", Astrophys. J. Lett., 547, L37-L40, (2001). 

[91] Edwards, R., and Bailes, M., "Recycled Pulsars Discovered at High Radio Frequency", As- 
trophys. J., 553, 801-808, (2001). 

[92] Edwards, R.T. , "Discovery of Eight Recycled Pulsars - The Swinburne Intermediate Latitude 
Pulsar Survey", in Kramer, M., Wex, N., and Wielebinski, R., eds., Pulsar Astronomy: 2000 
and Beyond (IAU Colloquium 177), Proceedings of the 177th Colloquium of the IAU, held 
at the Max-Planck-Institut fur Radioastronomie, Bonn, Germany, 30 August - 3 September 
1999, vol. 202 of ASP Conference Series, 33-34, (Astronomical Society of the Pacific, San 
Francisco, U.S.A., 2000). 



59 



[93] Faulkner, A.J., Kramer, M., Lyne, A.G., Manchester, R.N., McLaughlin, M.A., Stairs, I.H., 
Hobbs, G.B., Possenti, A., Lorimer, D.R., D'Amico, N, Camilo, F., and Burgay, M., "PSR 
J1756-2251: A New Relativistic Double Neutron Star System", Astrophys. J. Lett., 618, 
L119-L122, (2005). 

[94] Faulkner, A.J., Stairs, I.H., Kramer, M., Lyne, A.G., Hobbs, G.B., Possenti, A., Lorimer, 
D.R., Manchester, R.N., McLaughlin, M.A., D'Amico, N., Camilo, F., and Burgay, M., "The 
Parkes Multibeam Pulsar Survey - V. Finding binary and millisecond pulsars", Mon. Not. 
R. Astron. Soc, 355, 147-158, (2004). 

[95] Flannery, B.P., and van den Heuvel, E.P.J. , "On the origin of the binary pulsar PSR 
1913+16", Astron. Astrophys., 39, 61-67, (1975). 

[96] Fomalont, E.B., Goss, W.M., Lyne, A.G., Manchester, R.N., and Justtanont, K., "Positions 
and Proper Motions of Pulsars", Mon. Not. R. Astron. Soc, 258, 497-510, (1992). 

[97] Foster, R.S., and Backer, D.C., "Constructing a Pulsar Timing Array", Astrophys. J., 361, 
300-308, (1990). 

[98] Foster, R.S., Wolszczan, A., and Camilo, F., "A new binary millisecond pulsar", Astrophys. 
J. Lett, 410, L91-L94, (1993). 

[99] Freirc, P.C., Camilo, F., Kramer, M., Lorimer, D.R., Lyne, A.G., Manchester, R.N., and 
D'Amico, N., "Further results from the timing of the millisecond pulsars in 47 Tucanae", 
Mon. Not. R. Astron. Soc., 340, 1359-1374, (2003). 

[100] Freire, P.C., Gupta, Y., Ransom, S.M., and Ishwara-Chandra, C.H., "Giant Metrewave Radio 
Telescope Discovery of a Millisecond Pulsar in a Very Eccentric Binary System" , Astrophys. 
J. Lett., 606, L53-L56, (2004). 

[101] Freire, P.C.C., Hessels, J.W.T., Nice, D.J., Ransom, S.M., Lorimer, D.R., and Stairs, I.H., 
"The Millisecond Pulsars in NGC 6760", Astrophys. J., 621, 959-965, (2005). 

[102] Freire, P.C.C., Kramer, M., and Lyne, A.G., "Determination of the orbital parameters of 
binary pulsars", Mon. Not. R. Astron. Soc., 322, 885, (2001). Related online version (cited 

on 19 March 2001): 

http : //arXiv . org/abs/astro-ph/0010463 

[103] Fruchter, A.S., Stinebring, D.R., and Taylor, J.H., "A millisecond pulsar in an eclipsing 
binary", Nature, 333, 237-239, (1988). 

[104] Gil, J., and Melikidze, G.I., "On the Giant Nano-subpulses in the Crab Pulsar", in Camilo, 
F., and Gaensler, B.M., eds., Young neutron stars and their environments, Proceedings of 
an IAU Symposium on 14-17 July 2003 in Sydney, Australia, vol. 218 of ASP Conference 
Series, 321, (Astronomical Society of the Pacific, San Francisco, U.S.A., 2004). 

[105] Gold, T., "Rotating neutron stars as the origin of the pulsating radio sources", Nature, 218, 
731-732, (1968). 

[106] Gold, T., "Rotating neutron stars and the nature of pulsars", Nature, 221, 25-27, (1969). 

[107] Gould, D.M., and Lyne, A.G., "Multifrequency polarimetry of 300 radio pulsars" , Mon. Not. 
R. Astron. Soc, 301, 235-260, (1998). 

[108] Gunn, J.E., and Ostriker, J. P., "On the nature of pulsars. III. Analysis of observations", 
Astrophys. J., 160, 979-1002, (1970). 



60 



[109] Han, J.L., and Manchester, R.N., "The shape of radio pulsar beams", Mon. Not. R. Astron. 
Soc, 320, L35-L40, (2001). Related online version (cited on 19 March 2001): 
http : //arXiv . org/abs/astro-ph/0010538 

[110] Hankins, T.H., "Microsecond Intensity Variation in the Radio Emission from CP 0950", 
Astrophys. J., 169, 487-494, (1971). 

[Ill] Hankins, T.H., Kern, J.S., Weatherall, J.C., and Eilek, J. A., "Nanosecond radio bursts from 
strong plasma turbulence in the Crab pulsar", Nature, 422, 141-143, (2003). 

[112] Hannover University, "GEO 600 home page", project homepage, (2005). URL (cited on 29 
January 2005): 

http : //www . geo600 . uni-hannover . de 

[113] Hansen, B., and Phinney, E.S., "The pulsar kick velocity distribution", Mon. Not. R. Astron. 
Soc, 291, 569-577, (1997). 

[114] Hansen, B., and Phinney, E.S., "Stellar forensics - II. Millisecond pulsar binaries", Mon. 
Not. R. Astron. Soc, 294, 569-581, (1998). 

[115] Harrison, P.A., Lyne, A.G., and Anderson, B., "New Determinations of the Proper Motions 
of 44 Pulsars", Mon. Not. R. Astron. Soc, 261, 113-124, (1993). 

[116] Harvard-Smithsonian Center for Astrophysics, "The NASA Astrophysics Data System", web 
interface to database, (2001). URL (cited on 19 March 2001): 
http : // adsabs . harvard . edu/ 

[117] Hellings, R.W., "Thoughts on Detecting a Gravitational Wave Background with Pulsar 
Data", in Backer, D.C., ed., Impact of Pulsar Timing on Relativity and Cosmology, Proceed- 
ings of the workshop, held at Berkeley, June 7-9, 1990, kl, (Center for Particle Astrophysics, 
Berkeley, U.S.A., 1990). 

[118] Hellings, R.W., and Downs, G.S., "Upper limits on the isotropic gravitational radiation 
background from pulsar timing analysis", Astrophys. J. Lett., 265, L39-L42, (1983). 

[119] Hewish, A., Bell, S.J., Pilkington, J.D.H., Scott, P.F., and Collins, R.A., "Observation of a 
rapidly pulsating radio source", Nature, 217, 709-713, (1968). 

[120] Hills, J.G., "The effects of Sudden Mass loss and a random kick velocity produced in a 
supernova explosion on the Dynamics of a binary star of an Arbitrary orbital Eccentricity" , 
Astrophys. J., 267, 322-333, (1983). 

[121] Hobbs, G., "Pulsars and Gravitational Wave Detection", (2005). URL (cited on 22 January 

2005): 

http : //arXiv . org/abs/astro-ph/0412153 

[122] Hobbs, G., Lorimer, D.R., Lyne, A.G., and Kramer, M., "A statistical study of 233 pulsar 
proper motions", Mon. Not. R. Astron. Soc, 360, 974-992, (2005). 

[123] Hobbs, G.B., Faulkner, A., Stairs, I.H., Camilo, F., Manchester, R.N., Lyne, A.G., Kramer, 
M., D'Amico, N., Kaspi, V.M., Possenti, A., McLaughlin, M.A., Lorimer, D.R., Burgay, M., 
Joshi, B.C., and Crawford, F., "The Parkes multibeam pulsar survey: IV. Discovery of 176 
pulsars and parameters for 248 previously known pulsars", Mon. Not. R. Astron. Soc, 352, 
1439-1472, (2004). 



61 



[124] Hobbs, G.B., Lyne, A.G., Kramer, M., Martin, C.E., and Jordan, C, "Long-Term Timing 
Observations of 374 Pulsars", Mon. Not. R. Astron. Soc, 353, 1311-1344, (2004). 

[125] Hogan, C.J., and Rees, M.J., "Gravitational interactions of cosmic strings", Nature, 311, 
109-114, (1984). 

[126] Hotan, A.W., "Latest timing results for PSR J0437-4715", personal communication. 

[127] Hotan, A.W., Bailes, M., and Ord, S.M., "Geodetic Precession in PSR J1141-6545", Astro- 
phys. J. , 624, 906-913, (2005) . Related onlin e version (cited on 22 January 2005): 
http : //arXiv . org/abs/astro-ph/0412152 

[128] Hulse, R.A., "The discovery of the binary pulsar", Rev. Mod. Phys., 66, 699-710, (1994). 

[129] Hulse, R.A., and Taylor, J.H., "Discovery of a pulsar in a binary system", Astrophys. J. 
Lett, 195, L51-L53, (1975). 

[130] Hunt, G.C., "The rate of change of period of the pulsars", Mon. Not. R. Astron. Soc, 153, 
119-131, (1971). 

[131] I NFN, "The Virgo Project" , p roject homepage, (2005). URL (cited on 29 January 2005): 
http : //www .virgo . inf n. it 

[132] International SKA Project Office, "The Square Kilometre Array", project homepage, (2005). 
URL (cited on 23 January 2005): 
http : //www. skatelescope Torgl 

[133] Jacoby, B.A., Bailes, M., van Kerkwijk, M.H., Ord, S., Hotan, A., Kulkarni, S.R., and 
Anderson, S.B., "PSR J1909-3744: A Binary Millisecond Pulsar with a Very Small Duty 
Cycle", Astrophys. J. Lett, 599, L99-L102, (2003). 

[134] Jaffe, A.H., and Backer, D.C., "Gravitational Waves Probe the Coalescence Rate of Massive 
Black Hole Binaries", Astrophys. J., 583, 616-631, (2003). 

[135] Jenet, F.A., Lommen, A.N., Larson, S.L., and Wen, L., "Constraining the properties fo 
supermassive black hole systems using pulsar timing: Application to 3C 66B", Astrophys. 
J., 606, 799-803, (2004). 

[136] Johnston, H.M., and Kulkarni, S.R., "On the detectability of pulsars in close binary systems", 
Astrophys. J., 368, 504-514, (1991). 

[137] Johnston, S., "Evidence for a deficit of pulsars in the inner Galaxy", Mon. Not. R. Astron. 
Soc, 268, 595-601, (1994). 

[138] Johnston, S., Lorimer, D.R., Harrison, P.A., Bailes, M., Lyne, A.G., Bell, J.F., Kaspi, V.M., 
Manchester, R.N., D'Amico, N., Nicastro, L., and Shengzhen, J., "Discovery of a very bright, 
nearby binary millisecond pulsar", Nature, 361, 613-615, (1993). 

[139] Johnston, S., Lyne, A.G., Manchester, R.N., Kniffen, D.A., D'Amico, N., Lim, J., and 
Ashworth, M., "A High Frequency Survey of the Southern Galactic Plane for Pulsars" , Mon. 
Not. R. Astron. Soc, 255, 401-411, (1992). 

[140] Johnston, S., Manchester, R.N., Lyne, A.G., Bailes, M., Kaspi, V.M., Qiao, G., and D'Amico, 
N., "PSR 1259-63: A binary radio pulsar with a Be star companion", Astrophys. J. Lett., 
387, L37-L41, (1992). 



62 



[141] Jones, A.W., and Lyne, A.G., "Timing Observations of the Binary Pulsar PSR 0655+64", 
Mon. Not. R. Astron. Soc, 232, 473, (1988). 

[142] Joshi, B.C., Lyne, A.G., Kramer, M., Lorimer, D.R., Jordan, C.A., Holloway, A.J., Ikin, T.S., 
and Stairs, I.H., "Coherent On-line Baseband Receiver for Astronomy", in Bailes, M., Nice, 
D.J., and Thorsett, S.E., eds., Radio Pulsars, Proceedings of meeting held at Mediterranean 
Agonomic Institute of Chania, Crete, Greece, 26-29 August 2002, vol. 302 of ASP Conference 
Series, 321-322, (Astronomical Society of the Pacific, San Fransisco, U.S.A., 2003). 

[143] Jouteux, S.j Ramachandran, R., Stappers, B.W., Jonker, P.G., and van der Klis, M., "Search- 
ing for pulsars in close circular binary systems", Astron. Astrophys., 384, 532-544, (2002). 

[144] Kaaret, P., Morgan, E., and Vanderspek, R., "Orbital parameters of HETE J1900. 1-2455" , 
Astron. Telegram, 2005(538), (June, 2005). URL (cited on 19 July 2005): 
http : //www . astronomerstelegram . org/?read=538 

[145] Kalogera, V., "Close Binaries with Two Compact Objects", in Kramer, M., Wex, N., and 
Wiclebinski, R., eds., Pulsar Astronomy: 2000 and Beyond (IAU Colloquium 177), Pro- 
ceedings of the 177th Colloquium of the IAU, held at the Max-Planck-Institut fur Radioas- 
tronomie, Bonn, Germany, 30 August - 3 September 1999, ASP Conference Series, 579-584, 
(Astronomical Society of the Pacific, San Francisco, U.S.A., 2000). 

[146] Kalogera, V., "The Strongly Relativistic Double Pulsar and LISA", lecture notes, North- 
western University (2005). URL (cited on 19 July 2005): 

http : // www .astro . northwestern . edu/Vicky/TALKS/LISA_0737.ppt Invited talk at the 
5th International LISA Symposium, ESTEC, Noordwijk, The Netherlands, July 12-15, 2004. 

[147] Kalogera, V., Kim, C, Lorimer, D.R., Burgay, M., D'Amico, N., Possenti, A., Manchester, 
R.N., Lyne, A.G., Joshi, B.C., McLaughlin, M.A., Kramer, M., Sarkissian, J.M., and Camilo, 
F., "The Cosmic Coalescence Rates for Double Neutron Star Binaries", Astrophys. J. Lett., 
601, L179-L182, (2004). 

[148] Kalogera, V., Kim, C, Lorimer, D.R., Burgay, M., D'Amico, N., Possenti, A., Manchester, 
R.N., Lyne, A.G., Joshi, B.C., McLaughlin, M.A., Kramer, M., Sarkissian, J.M., and Camilo, 
F., "Erratum: "The Cosmic Coalescence Rates for Double Neutron Star Binaries"", Astro- 
phys. J. Lett., 614, L137-L138, (2004). 

[149] Kalogera, V., Kim, C, Lorimer, D.R., Ihm, M., and Belczynski, C, "The Galactic formation 
rate of eccentric neutron star white dwarf binaries", in Rasio, F.A., and Stairs, I.H., eds., 
Binary Radio Pulsars, Proceedings of meeting held at the Aspen Center for Physics, USA, 
12-16 January 2004, vol. 328 of ASP Conference Series, 261-267, (Astronomical Society of 
the Pacific, San Francisco, U.S.A., 2005). 

[150] Kalogera, V., and Lorimer, D. R., "An Upper Limit on the Coalescence Rate of Double 
Neutron-Star Binaries in the Galaxy", Astrophys. J., 530, 890-895, (2000). 

[151] Kalogera, V., Narayan, R., Spergel, D.N., and Taylor, J.H., "The Coalescence Rate of Double 
Neutron Star Systems", Astrophys. J., 556, 340-356, (2001). Related online version (cited 

on 19 March 2001): 

http : //arXiv . org/abs/astro-ph/0012038 

[152] Kaspi, V.M., Applications of Pulsar Timing, Ph.D. Thesis, (Princeton University, Princeton, 
U.S.A., 1994). 



63 



[153] Kaspi, V.M., Johnston, S., Bell, J.F., Manchester, R.N., Bailes, M., Bessell, M., Lync, A.G., 
and D'Amico, N., "A massive radio pulsar binary in the Small Magellanic Cloud", Astrophys. 
J. Lett, 423, L43-L45, (1994). 

[154] Kaspi, V.M., Lyne, A.G., Manchester, R.N., Crawford, F., Camilo, F., Bell, J.F., D'Amico, 
N., Stairs, I.H., McKay, N.P.F., Morris, D.J., and Possenti, A., "Discovery of a young radio 
pulsar in a relativistic binary orbit", Astrophys. J., 534, 321-327, (2000). 

[155] Kaspi, V.M., Taylor, J.H., and Ryba, M.F., "High-Precision Timing of Millisecond Pulsars. 
III. Long-Term Monitoring of PSRs B1855+09 and B1937+21" , Astrophys. J., 428, 713-728, 
(1994). 

[156] Kim, C, Kalogera, V., and Lorimer, D.R., "The Probability Distribution of Binary Pulsar 
Coalescence Rates. I. Double Neutron Star Systems in the Galactic Field", Astrophys. J., 
584, 985-995, (2003). 

[157] Kim, C, Kalogera, V., Lorimer, D.R., and White, T., "The Probability Distribution of 
Binary Pulsar Coalescence Rates. II. Neutron Star- White Dwarf Binaries", Astrophys. J., 
616, 1109-1117, (2004). 

[158] Kluzniak, W., Ruderman, M., Shaham, J., and Tavani, M., "Nature and evolution of the 
eclipsing millisecond binary pulsar PSR 1957+20", Nature, 334, 225-227, (1988). 

[159] Kolonko, M., Gil, J., and Maciesiak, K., "On the pulse-width statistics in radio pulsars", 
Astron. Astrophys., 428, 943-951, (2005). 

[160] Kopcikin, S.M., "Binary pulsars as detectors of ultralow-frequency gravitational waves", 
Phys. Rev. D, 56, 4455-4469, (1997). 

[161] Kramer, M., "Latest timing results for PSR J0737-3039", personal communication. 

[162] Kramer, M., Backer, D.C., Cordes, J.M., Lazio, T.J.W., Stappcrs, B.W., and Johnston, S., 
"Strong-field tests of gravity using pulsars and black holes" , New Astron. Rev., 48, 993-1002, 
(2004). 

[163] Kramer, M., Bell, J.F., Manchester, R.N., Lyne, A.G., Camilo, F., Stairs, I.H., D'Amico, N., 
Kaspi, V.M., Hobbs, G., Morris, D.J., Crawford, F., Possenti, A., Joshi, B.C., McLaughlin, 
M.A., Lorimer, D.R., and Faulkner, A.J., "The Parkes Multibeam Pulsar Survey - III. Young 
pulsars and the discovery and timing of 200 pulsars", Mon. Not. R. Astron. Soc, 342, 1299- 
1324, (2003). 

[164] Kramer, M., Xilouris, K.M., Camilo, F., Nice, D.J., Backer, D.C., Lorimer, D.R., 
Doroshenko, O.V., Lange, C, and Sallmen, S., "Profile Instabilities of the Millisecond Pulsar 
J1022+1001", Astrophys. J., 520, 324-334, (1999). 

[165] Kramer, M., Xilouris, K.M., Lorimer, D.R., Doroshenko, O.V., Jessner, A., Wielebinski, 
R., Wolszczan, A., and Camilo, F., "The characteristics of millisecond pulsar emission: I. 
Spectra, pulse shapes and the beaming fraction", Astrophys. J., 501, 270-285, (1998). 

[166] Kulkarni, S.R., Anderson, S.B., Prince, T.A., and Wolszczan, A., "Old pulsars in the low- 
density globular clusters M13 and M53", Nature, 349, 47-49, (1991). 

[167] Lange, C, Camilo, F., Wex, N., Kramer, M., Backer, D.C., Lyne, A.G., and Doroshenko, 
O.V., "Precision timing measurements of PSR J1012+5307", Mon. Not. R. Astron. Soc., 
326, 274-282, (2001). 



64 



[168] Lattimer, J.M., and Prakash, M., "The Physics of Neutron Stars", Science, 304, 536-542, 
(2004). 



[169] Lawson, K.D., Mayer, C.J., Osborne, J.L., and Parkinson, M.L., "Variations in the Spectral 
Index of the Galactic Radio Continuum Emission in the Northern Hemisphere" , Mon. Not. 
R. Astron. Soc, 225, 307-327, (1987). 

[170] Levinson, A., and Eichler, D., "Can neutron stars ablate their companions?", Astrophys. J., 
379, 359-365, (1991). 

[171] Lewandowski, W., Wolszczan, A., Feiler, C, Konacki, M., and Soltysihski, T., "Arecibo tim- 
ing and single-pulse observations of eighteen pulsars", Astrophys. J., 600, 905-913, (2004). 

[172] Li, X., "Formation of Mildly Recycled Low-Mass Binary Pulsars from Intermediate-Mass 
X-Ray Binaries", Astrophys. J., 564, 930-934, (2002). 

[173] Lipunov, V.M., Postnov, K.A., and Prokhorov, M.E., "The Scenario Machine: Binary Star 
Population Synthesis", Astrophys. Space Phys. Rev., 9, 1-178, (1996). Related online version 

(cited on 24 October 2005): 

http : //xray . sai .msu.ru/~mystery/articles/review/ 

[174] Lohmer, O., Kramer, M., Driebe, T., Jessner, A., Mitra, D., and Lyne, A.G., "The parallax, 
mass and age of the PSR J2145-0750 binary system", Astron. Astrophys., 426, 631-640, 
(2004). 

[175] Lommen, A.N., "New Limits on Gravitational Radiation using Pulsars", in Becker, W., 
Lesch, H., and Truemper, J., eds., Neutron Stars, Pulsars and Supernova Remnants, Pro- 
ceedings of 270. WE-Heraeus Seminar 21-25 January 2002, Physikzcntrum Bad Honnef, Ger- 
many, vol. 278 of MPE-Report, 114-125, (Max Planck Institute for Extraterrestrial Physics, 
Garching, Germany, 2002). Related online version (cited on 24 January 2005): 
ftp : //ftp . xray .mpe .mpg . de/people/web/MPE_Report_278 . pdf 

[176] Lommen, A.N., and Backer, D.C., "Using pulsars to detect massive black hole binaries via 
gravitational radiation: Sagittarius A* and nearby galaxies", Astrophys. J., 562, 297-302, 
(2001). 

[177] Lommen, A.N., Zepka, A.F., Backer, D.C., McLaughlin, M.A., Cordes, J.M., Arzoumanian, 
Z., and Xilouris, K.M., "New Pulsars from an Arecibo Drift Scan Search", Astrophys. J., 
545, 1007-1014, (2001). 

[178] Lorimer, D.R., "Pulsar statistics - II. The local low-mass binary pulsar population", Mon. 
Not. R. Astron. Soc, 274, 300-304, (1995). 

[179] Lorimer, D.R., "Binary and Millisecond Pulsars", Living Rev. Relativity, 1, lrr-1998-10, 

(1998). URL (cited on 19 March 2001): 

http : //www . livingreviews . org/lrr- 1998- 10 

[180] Lorimer, D.R., "Binary and Millisecond Pulsars", Living Rev. Relativity, 4, lrr-2001-5, 
(2001). URL (cited on 22 January 2005): 
http : //www . livingreviews . org/lrr- 2001- 5 

[181] Lorimer, D.R., "The Galactic Population and Birth Rate of Radio Pulsars", in Camilo, F., 
and Gaensler, B.M., eds., Young neutron stars and their environments, Proceedings of an 
IAU Symposium on 14-17 July 2003 in Sydney, Australia, vol. 218 of ASP Conference Series, 
105, (Astronomical Society of the Pacific, San Francisco, U.S.A., 2004). 



65 



[182] Lorimer, D.R., Bailes, M., Dewey, R.J., and Harrison, P.A., "Pulsar statistics: the birthrate 
and initial spin periods of radio pulsars", Mon. Not. R. Astron. Soc, 263, 403-415, (1993). 



[183] Lorimer, D.R., Camilo, F., Freire, P.C.C., Kramer, M., Lyne, A.G., Manchester, R.N., and 
D'Amico, N., "Millisecond radio pulsars in 47 Tucanae", in Bailes, M., Nice, D.J., and 
Thorsett, S.E., eds., Radio Pulsars, Proceedings of meeting held at Mediterranean Agonomic 
Institute of Chania, Crete, Greece, 26-29 August 2002, vol. 302 of ASP Conference Series, 
363-366, (Astronomical Society of the Pacific, San Fransisco, U.S.A., 2003). 

[184] Lorimer, D.R., and Freire, P.C.C., "Long-period binaries and the strong equivalence prin- 
ciple", in Rasio, FA., and Stairs, I.H., eds., Binary Radio Pulsars, Proceedings of meeting 
held at the Aspen Center for Physics, USA, 12-16 January 2004, vol. 328 of ASP Conference 
Series, 19-24, (Astronomical Society of the Pacific, San Francisco, U.S.A., 2005). 

[185] Lorimer, D.R., and Kramer, M., Handbook of Pulsar Astronomy, vol. 4 of Cambridge Observ- 
ing Handbooks for Research Astronomers, (Cambridge University Press, Cambridge, U.K.; 
New York, U.S.A, 2005), 1st edition. 

[186] Lorimer, D.R., and Kramer, M., Handbook of Pulsar Astronomy, vol. 4 of Cambridge Observ- 
ing Handbooks for Research Astronomers, (Cambridge University Press, Cambridge, U.K.; 
New York, U.S.A., 2005). Related online version (cited on 18 January 2005): 
http : //www . jb .mail. ac .uk/pulsarhandbook 

[187] Lorimer, D.R., Lyne, A.G., Bailes, M., Manchester, R.N., D'Amico, N., Stappers, B.W., 
Johnston, S., and Camilo, F., "Discovery of four binary millisecond pulsars", Mon. Not. R. 
Astron. Soc, 283, 1383-1387, (1996). 

[188] Lorimer, D.R., McLaughlin, M.A., Arzoumanian, Z., Xilouris, K.M., Cordes, J.M., Lommen, 
A.N., Fruchter, A.S., Chandler, A.M., and Backer, D.C., "PSR J0609+2130: a disrupted 
binary pulsar?", Mon. Not. R. Astron. Soc, 347, L21-L25, (2004). 

[189] Lorimer, D.R., Nicastro, L., Lyne, A.G., Bailes, M., Manchester, R.N., Johnston, S., Bell, 
J.F., DAmico, N., and Harrison, PA., "Four new millisecond pulsars in the Galactic disk", 
Astrophys. J., 439, 933-938, (1995). 

[190] Lorimer, D.R., Xilouris, KM., Fruchter, A.S., Stairs, I.H., Camilo, F., McLaughlin, M.A., 
Vazquez, A.M., Eder, J.E., Roberts, M.S.E., Hessels, J.W.T., and Ransom, S.M., "Discovery 
of 10 pulsars in an Arecibo drift-scan survey", Mon. Not. R. Astron. Soc, 359, 1524-1530, 
(2005). Related online version (cited on 19 July 2005): 
http : //arXiv . org/abs/astro-ph/0504019 

[191] Lorimer, D.R., Yates, J. A., Lyne, A.G., and Gould, DM., "Multifrequency Flux density 
measurements of 280 pulsars", Mon. Not. R. Astron. Soc, 273, 411-421, (1995). 

[192] Lundgren, S.C., Zepka, A.F., and Cordes, J.M., "A millisecond pulsar in a six hour orbit: 
PSR J0751+1807", Astrophys. J., 453, 419-423, (1995). 

[193] Lyne, A.G., "Binary Pulsar Discoveries in the Parkes Multibeam Survey", in Rasio, FA., 
and Stairs, I.H., eds., Binary Radio Pulsars, Proceedings of meeting held at the Aspen Center 
for Physics, USA, 12 Janaury - 16 January 2004, vol. 328 of ASP Conference Series, 37-41, 
(Astronomical Society of the Pacific, San Francisco, U.S.A., 2005). 

[194] Lyne, A.G., Anderson, B., and Salter, M.J., "The Proper Motions of 26 Pulsars", Mon. Not. 
R. Astron. Soc, 201, 503-520, (1982). 



66 



[195] Lyne, A.G., and Bailes, M., "The mass of the PSR 2303+46 system", Mon. Not. R. Astron. 
Soc, 246, 15P-17P, (1990). 

[196] Lyne, A.G., Biggs, J.D., Harrison, P.A., and Bailes, M., "A long-period globular-cluster 
pulsar in an eclipsing binary system", Nature, 361, 47-49, (1993). 

[197] Lyne, A.G., Brinklow, A., Middleditch, J., Kulkarni, S.R., Backer, D.C., and Clifton, T.R., 
"The discovery of a millisecond pulsar in the globular cluster M28", Nature, 328, 399-401, 
(1987). 

[198] Lyne, A.G., Burgay, M., Kramer, M., Possenti, A., Manchester, R.N., Camilo, F., McLaugh- 
lin, M., Lorimer, D.R., Joshi, B.C., Reynolds, J.E., and Freire, P.C.C., "A Double-Pulsar 
System: A Rare Laboratory for Relativistic Gravity and Plasma Physics", Science, 303, 
1153-1157, (2004). 

[199] Lyne, A.G., Camilo, F., Manchester, R.N., Bell, J.F., Kaspi, V.M., D'Amico, N., McKay, 
N.P.F., Crawford, F., Morris, D.J., Sheppard, D.C., and Stairs, I.H., "The Parkes multibeam 
survey: PSR J1811-1736 - a pulsar in a highly eccentric binary system", Mon. Not. R. 
Astron. Soc, 312, 698-702, (2000). 

[200] Lyne, A.G., and Graham-Smith, F., Pulsar Astronomy, Cambridge Astrophysics Series, 
(Cambridge University Press, Cambridge, U.K.; New York, U.S. A, 2006), 3rd edition, in 
press. 

[201] Lyne, A.G., and Lorimer, D.R., "High birth velocities of radio pulsars", Nature, 369, 127- 
129, (1994). 

[202] Lyne, A.G., and Manchester, R.N., "The shape of pulsar radio beams", Mon. Not. R. Astron. 
Soc, 234, 477-508, (1988). 

[203] Lyne, A.G., Manchester, R.N., and D'Amico, N., "PSR B1745-20 and young pulsars in 
globular clusters", Astrophys. J. Lett., 460, L41-L44, (1996). 

[204] Lyne, A.G., Manchester, R.N., D'Amico, N., Staveley-Smith, L., Johnston, S., Lim, J., 
Fruchter, A.S., Goss, W.M., and Frail, D.A., "An eclipsing millisecond pulsar in the globular 
cluster Terzan 5", Nature, 347, 650-652, (1990). 

[205] Lyne, A.G., Manchester, R.N., Lorimer, D.R., Bailes, M., D'Amico, N., Tauris, T.M., John- 
ston, S., Bell, J.F., and Nicastro, L., "The Parkes Southern Pulsar Survey - II. Final Results 
and Population Analysis", Mon. Not. R. Astron. Soc, 295, 743-755, (1998). 

[206] Lyne, A.G., Mankelow, S.H., Bell, J.F., and Manchester, R.N., "Radio pulsars in Terzan 5", 
Mon. Not. R. Astron. Soc, 316, 491-493, (2000). 

[207] Lyne, A.G., and McKenna, J., "PSR 1820 - 11: a binary pulsar in a wide and highly eccentric 
orbit", Nature, 340, 367-369, (1989). 

[208] Manchester, R.N., Hobbs, G.B., Teoh, A., and Hobbs, M., "The ATNF Pulsar Catalogue", 
(2004). URL (cited on 22 January 2005): 
http : //arXiv . org/abs/astro-ph/0412641 

[209] Manchester, R.N., and Lyne, A.G., "Pulsar Interpulses - two poles or one?", Mon. Not. R. 
Astron. Soc, 181, 761-767, (1977). 



67 



[210] Manchester, R.N., Lyne, A.G., Camilo, F., Bell, J.F., Kaspi, V.M., D'Amico, N., McKay, 
N.P.F., Crawford, F., Stairs, I.H., Possenti, A., Morris, D.J., and Sheppard, D.C., "The 
Parkes multi-beam pulsar survey - I. Observing and data analysis systems, discovery and 
timing of 100 pulsars", Mon. Not. R. Astron. Soc., 328, 17-35, (2001). 

[211] Manchester, R.N., Lyne, A.G., D'Amico, N., Bailes, M., Johnston, S., Lorimer, D.R., Har- 
rison, P.A., Nicastro, L., and Bell, J.F., "The Parkes Southern Pulsar Survey I. Observing 
and data analysis systems and initial results", Mon. Not. R. Astron. Soc, 279, 1235-1250, 
(1996). 

[212] Manchester, R.N., Newton, L.M., Cooke, D.J., and Lyne, A.G., "Detection of a pulsar in a 
long-period binary system", Astrophys. J. Lett, 236, L25-L27, (1980). 

[213] Manchester, R.N., and Taylor, J.H., "Period irregularities in pulsars", Astrophys. J. Lett., 
191, L63-L65, (1974). 

[214] Matsakis, D.N., Taylor, J.H., Eubanks, T.M., and Marshall, T., "A statistic for describing 
pulsar and clock stabilities", Astron. Astrophys., 326, 924-928, (1997). 

[215] Max Planck Institute for Radio Astronomy, "MPIfR Pulsar Group", project homepage, 

(1999). URL (cited on 20 March 2001): 

http : //www . mpif r-bonn . mpg . de/div/pulsar/ 

[216] McLaughlin, M.A., Kramer, M., Lyne, A.G., Lorimer, D.R., Stairs, I.H., Possenti, A., Manch- 
ester, R.N., Freire, P.C.C., Joshi, B.C., Burgay, M., Camilo, F., and DAmico, N., "The 
double pulsar system J0737-3039: Modulation of the radio emission from B by the radiation 
from A", Astrophys. J. Lett, 613, L57-L60, (2004). 

[217] McLaughlin, M.A., Lorimer, D.R., Champion, D.J., Arzoumanian, Z., Backer, D.C., Cordes, 
J.M., Fruchter, A.S., Lommen, A.N., and Xilouris, K.M., "New Binary and Millisecond 
Pulsars from Arecibo Drift-Scan Searches", in Rasio, F.A., and Stairs, I.H., eds., Binary 
Radio Pulsars, Proceedings of meeting held at the Aspen Center for Physics, USA, 12 Janaury 
- 16 January 2004, vol. 328 of ASP Conference Series, 43-49, (Astronomical Society of the 
Pacific, San Francisco, U.S.A., 2005). 

[218] Michel, F.C., Theory of Neutron Star Magnetospheres, Theoretical Astrophysics, (University 
of Chicago Press, Chicago, U.S.A., 1991). 

[219] Middleditch, J., and Kristian, J., "A search for young, luminous optical pulsars in extra- 
galactic supernova remnants", Astrophys. J., 279, 157-161, (1984). 

[220] Morgan, E., Kaaret, P., and Vanderspek, R., "HETE J1900. 1-2455 is a millisecond pulsar", 
Astron. Telegram , 2005(5 23), (June, 20 05). URL (cited on 19 July 2005): 
http : //www . astronomerstelegram . org/?read=523 

[221] Morris, D.J., Hobbs, G., Lyne, A.G., Stairs, I.H., Camilo, F., Manchester, R.N., Possenti, 
A., Bell, J.F., Kaspi, V.M., DAmico, N., McKay, N.P.F., Crawford, F., and Kramer, M., 
"The Parkes Multibeam Pulsar Survey - II. Discovery and timing of 120 pulsars" , Mon. Not. 
R. Astron. Soc., 335, 275-290, (2002). 

[222] Narayan, R., "The birthrate and initial spin period of single radio pulsars", Astrophys. J., 
319, 162-179, (1987). 

[223] Narayan, R., and Vivekanand, M., "Evidence for evolving elongated Pulsar Beams", Astron. 
Astrophys., 122, 45-53, (1983). 



68 



[224] NASA/ESA, "Laser Interferometer Space Antenna", project homepage, (2005). URL (cited 
on 24 January 2005): 
http : //lisa . jpl . nasa . gov 

[225] National Astronomical Observatory, "TAMA project" , project homepage, (2005). URL (cited 

on 29 January 2005): 

http : / /t amago . mtk . n ao . ac . jp| 

[226] National Astronomy and Ionosphere Center, "Pulsars at the Arecibo Observatory" , project 
homepag e, (1999). URL (c ited on 7 December 2000): 
http : / / www . naic . edu/~pulsar 

[227] National Astronomy and Ionosphere Center, "Arecibo L-band feed array" , project homepage, 
(2005). URL (cited on 23 January 2005): 
http : //alf a Tnaic . edu| 

[228] National Astronomy and Ionosphere Center, "The Arecibo Observatory" , project homepage, 
(2005). URL (cited on 29 January 2005): 
http : // www . naic . edu 

[229] National Astronomy and Ionosphere Center, "Preliminary PALFA survey", project home- 
page, (2005). URL (cited on 23 January 2005): 
http : //www. naic . edu/~palf a 

[230] National Astronomy and Ionospheric Center, "Pulsars in globular clusters" , web interface to 
database, (2005). URL (cited on 18 January 2005): 
http : //www . naic . edu/~pf reire/GCpsr .html 

[231] National Centre for Radio Astonomy, "The Giant Metre Wave Radio Telescope", project 
homepage, (2005). URL (cited on 23 January 2005): 
|http : // www . gmr t .ncra. tif r .res . in| 

[232] National Radio Astronomy Observatory, "The Green Bank Telescope", project homepage, 
(2005). URL (cited on 23 J anuary 2005): 
http : // www . gb . nrao . edu 

[233] National Radio Astronomy Observatory, "Scott Ransom's home page" , personal homepage, 
(2005). URL (cited on 7 November 2005): 
http : / / www . cv . nrao . edu/~sransom 

[234] Navarro, J., Anderson, S.B., and Freire, P.C.C., "The Arecibo 430 MHz Intermediate Galactic 
Latitude Survey: Discovery of Nine Radio Pulsars", Astrophys. J., 594, 943-951, (2003). 

[235] Navarro, J., de Bruyn, A.G., Frail, D.A., Kulkarni, S.R., and Lyne, A.G., "A very luminous 
binary millisecond pulsar", Astrophys. J. Lett., 455, L55-L58, (1995). 

[236] Nicastro, L., Lyne, A.G., Lorimer, D.R., Harrison, P.A., Bailes, M., and Skidmore, B.D., 
"PSR J1012+5307: a 5.26-ms pulsar in a 14.5-h binary system", Mon. Not. R. Astron. Soc, 
273, L68-L70, (1995). 

[237] Nice, D.J., Sayer, R.W., and Taylor, J.H., "PSR J 1518+4904: A Mildly Relativistic Binary 
Pulsar System", Astrophys. J. Lett., 466, L87-L90, (1996). 



69 



[238] Nice, D.J., Sayer, R.W., and Taylor, J.H., "Timing Observations of the J1518+4904 Double 
Neutron Star System", in Arzoumanian, Z., van der Hooft, F., and van den Heuvel, E.P.J., 
eds., Pulsar Timing, General Relativity, and the Internal Structure of Neutron Stars, Pro- 
ceedings of the colloquium, Amsterdam, 24-27 September 1996, 79-83, (Royal Netherlands 
Academy of Arts and Sciences, Amsterdam, Netherlands, 1999). 

[239] Nice, D.J., Splaver, E.M., and Stairs, I.H., "On the Mass and Inclination of the PSR 
J2019+2425 Binary System", Astrophys. J., 549, 516-521, (2001). 

[240] Nice, D.J., and Taylor, J.H., "PSRs J2019+2425 and J2322+2057 and the Proper Motions 
of Millisecond Pulsars", Astrophys. J., 441, 429-435, (1995). 

[241] Ord, S.M., Bailes, M., and van Straten, W., "The Scintillation Velocity of the Relativistic 
Binary Pulsar PSR J1141-6545", Astrophys. J. Lett., 574, L75-L78, (2002). 

[242] Pacini, F., "Rotating neutron stars, pulsars, and supernova remnants" , Nature, 219, 145-146, 
(1968). 

[243] Peebles, P.J.E., Principles of Physical Cosmology, Princeton Series in Physics, (Princeton 
University Press, Princeton, U.S.A., 1993). 

[244] Penn State University, "Pulsar Planets" , personal homepage, (2001). URL (cited on 19 March 

2001): 

http : //www . astro . psu . edu/users/alex/pulsar_planets . htm 

[245] Peters, P.C., "Gravitational Radiation and the Motion of Two Point Masses", Phys. Rev., 
136, B1224-B1232, (1964). 

[246] Peters, P.C., and Mathews, J., "Gravitational Radiation from Point Masses in a Keplerian 
Orbit", Phys. Rev., 131, 435-440, (1963). 

[247] Pfahl, E., and Loeb, A., "Probing the Spacetime Around Sgr A* with Radio Pulsars", 
Astrophys. J., 615, 253-258, (2004). 

[248] Phinney, E.S., "The rate of neutron star binary mergers in the universe: Minimal predictions 
for gravity wave detectors", Astrophys. J. Lett., 380, L17-L21, (1991). 

[249] Phinney, E.S., "Pulsars as Probes of Newtonian Dynamical Systems", Philos. Trans. R. Soc. 
London, 341, 39-75, (1992). 

[250] Phinney, E.S., and Blandford, R.D., "Analysis of the Pulsar P-P distribution", Mon. Not. 
R. Astron. Soc, 194, 137-148, (1981). 

[251] Phinney, E.S., and Kulkarni, S.R., "Binary and millisecond pulsars", Annu. Rev. Astron. 
Astrophys., 32, 591-639, (1994). 

[252] Phinney, E.S., and Sigurdsson, S., "Ejection of pulsars and binaries to the outskirts of glob- 
ular clusters", Nature, 349, 220-223, (1991). 

[253] Phinney, E.S., and Verbunt, F., "Binary pulsars before spin-up and PSR 1820 - 11", Mon. 
Not. R. Astron. Soc, 248, 21P-23P, (1991). 

[254] Portegies Zwart, S.F., and Verbunt, F., "Population synthesis of high-mass binaries" , Astron. 
Astrophys., 309, 179, (1996). 



70 



[255] Possenti, A., D'Amico, N., Conongiu, A., Manchester, R.N., Sarkissian, J.M., Camilo, F., 
and Lyne, A.G., "A dozen new pulsars from the Parkes globular cluster search", in Rasio, 
F.A., and Stairs, I.H., eds., Binary Radio Pulsars, Proceedings of meeting held at the Aspen 
Center for Physics, USA, 12-16 January 2004, vol. 328 of ASP Conference Series, 189-194, 
(Astronomical Society of the Pacific, San Francisco, U.S.A., 2005). 

[256] Possenti, A., D'Amico, N., Manchester, R.N., Camilo, F., Lyne, A.G., Sarkissian, J., and 
Corongiu, A., "Three Binary Millisecond Pulsars in NGC 6266", Astrophys. J., 599, 475-484, 
(2003). 

[257] Prince, T.A., Anderson, S.B., Kulkarni, S.R., and Wolszczan, A., "Timing observations of 
the 8 hour binary pulsar 2127+11C in the globular cluster M15", Astrophys. J. Lett., 374, 
L41-L44, (1991). 

[258] Princeton University, "Princeton Pulsar Group", project homepage, (1995). URL (cited on 

20 March 2001): 

http : //pulsar . princeton . edu/ 

[259] Pulsar Group of the Australia Telescope National Facility (ATNF), "Parkes Multibeam Pul- 
sar Survey", project homepage, (2000). URL (cited on 2 December 2000): 
http : //www. atnf . csiro . au/research/pulsar/pmsurv/ 

[260] Rajagopal, M., and Romani, R.W., "Ultra-low-frequency gravitatational radiation from mas- 
sive black hole binaries", Astrophys. J., 446, 543-549, (1995). 

[261] Rankin, J.M., "Toward an empirical theory of pulsar emission. I. Morphological taxonomy", 
Astrophys. J., 274, 333-358, (1983). 

[262] Ransom, S.M., New search techniques for binary pulsars, Ph.D. Thesis, (Harvard University, 
Harvard, U.S.A., 2001). 

[263] Ransom, S.M., Cordes, J.M., and Eikenberry, S.S., "A New Search Technique for Short 
Orbital Period Binary Pulsars", Astrophys. J. Lett, 589, 911-920, (2003). 

[264] Ransom, S.M., Greenhill, L., Herrnstein, J.R., Manchester, R.N., Camilo, F., Eikenberry, 
S.S., and Lyne, A.G., "A Binary Millisecond Pulsar in Globular Cluster NGC 6544", Astro- 
phys. J. Lett, 546, L25-L28, (2001). 

[265] Ransom, S.M., Hessels, J.H., Stairs, I.H., Kaspi, V.M., and Backer, D.C., "A 20-cm survey 
for pulsars in globular clusters using the GBT and Arecibo" , in Rasio, F.A., and Stairs, I.H., 
eds., Binary Radio Pulsars, Proceedings of meeting held at the Aspen Center for Physics, 
USA, 12-16 January 2004, vol. 328 of ASP Conference Series, 199-205, (Astronomical Society 
of the Pacific, San Francisco, U.S.A., 2005). 

[266] Ransom, S.M., Hessels, J.W.T., Stairs, I.H., Freire, P.C.C., Kaspi, V.M., Camilo, F., and Ka- 
plan, D.L., "Twenty-One Millisecond Pulsars in Terzan 5 Using the Green Bank Telescope", 
Science, 307(5711), 892-896, (2005). 

[267] Ransom, S.M., Stairs, I.H., Backer, D.C., Greenhill, L.J., Bassa, C.G., Hessels, J.W.T., and 
Kaspi, V.M., "GBT Discovery of Two Binary Millisecond Pulsars in the Globular Cluster 
M30", Astrophys. J., 604, 328, (2004). 

[268] Ratnatunga, K.U., and van den Bergh, S., "The rate of stellar collapses in the galaxy", 
Astrophys. J., 343, 713-717, (1989). 



71 



[269] Ray, P.S., Thorsett, S.E., Jenet, F.A., van Kerkwijk, M.H., Kulkarni, S.R., Prince, T.A., 
Sandhu, J.S., and Nice, D.J., "A Survey for Millisecond Pulsars", Astrophys. J., 470, 1103- 
1110, (1996). 

[270] Richards, D.W., and Cornelia, J.M., "The period of pulsar NP 0532", Nature, 222, 551-552, 
(1969). 

[271] Romani, R.W., "Timing a millisecond pulsar array", in Ogelman, H., and van den Heuvel, 
E.P.J. , cds., Timing Neutron Stars, Proceedings of the NATO Advanced Study Institute on 
Timing Neutron Stars, gesme, Izmir, Turkey, 4-15 April 1988, vol. 262 of NATO ASI Series, 
Series C, 113-117, (Kluwer, Dordrecht, Netherlands; Boston, U.S.A., 1989). 

[272] Romani, R.W., "Populations of low-mass black hole binaries", Astrophys. J., 399, 621-626, 
(1992). 

[273] Romani, R.W., and Taylor, J.H., "An Upper Limit on the Stochastic Background of Ultralow- 
Frequency Gravitational Waves", Astrophys. J. Lett., 265, L35-L37, (1983). 

[274] Ryba, M.F., and Taylor, J.H., "High Precision Timing of Millisecond Pulsars. I. Astrometry 
and Masses of the PSR 1855+09 System", Astrophys. J., 371, 739-748, (1991). 

[275] Ryba, M.F., and Taylor, J.H., "High Precision Timing of Millisecond Pulsars. II. Astrometry, 
Orbital Evolution, and Eclipses of PSR 1957+20", Astrophys. J., 380, 557-563, (1991). 

[276] Sandhu, J.S., Bailes, M., Manchester, R.N., Navarro, J., Kulkarni, S.R., and Anderson, S.B., 
"The Proper Motion and Parallax of PSR J0437-4715", Astrophys. J. Lett, 478, L95-L98, 
(1997). 

[277] Sazhin, M.V., "Opportunities for detecting ultralong gravitational waves" , Sov. Astron., 22, 
36-38, (1978). 

[278] Scheuer, P.A.G., "Amplitude variations of pulsed radio sources", Nature, 218, 920-922, 
(1968). 

[279] Segelstein, D.J., Rawley, L.A., Stinebring, D.R., Fruchter, A.S., and Taylor, J.H., "New 
millisecond pulsar in a binary system", Nature, 322, 714-717, (1986). 

[280] Seiradakis, J.H., and Wielebinski, R., "Morphology and characteristics of radio pulsars", 
Astron. Astrophys. Rev., 12, 239-271, (2004). 

[281] Shemar, S.L., and Lyne, A.G., "Observations of pulsar glitches", Mon. Not. R. Astron. Soc, 
282, 677-690, (1996). 

[282] Shklovskii, I.S., "Possible causes of the secular increase in pulsar periods", Sov. Astron., 13, 
562-565, (1970). 

[283] Sigurdsson, S., and Phinncy, E.S., "Dynamics and interactions of binaries and neutron stars 
in globular clusters", Astrophys. J. Suppl. Ser., 99, 609-635, (1995). 

[284] Sigurdsson, S., Richer, H.B., Hansen, B.M., Stairs, I.H., and Thorsett, S.E., "A Young White 
Dwarf Companion to Pulsar B1620-26: Evidence for Early Planet Formation", Science, 301, 
193-196, (2003). 

[285] Smarr, L.L., and Blandford, R.D., "The binary pulsar: Physical processes, possible compan- 
ions and evolutionary histories", Astrophys. J., 207, 574-588, (1976). 



72 



[286] Smith, F.G., "The radio emission from pulsars", Rep. Prog. Phys., 66, 173-238, (2003). 



[287] Splaver, E.M., Nice, D.J., Stairs, I.H., Lommen, A.N., and Backer, D.C., "Masses, Parallax, 
and Relativistic Timing of the PSR J1713+0747 Binary System", Astrophys. J., 620, 405- 
415, (2005). 

[288] Stairs, I. H., "Overview of Pulsar Tests of General Relativity" , in Rasio, F.A., and Stairs, I.H., 
eds., Binary Radio Pulsars, Proceedings of meeting held at the Aspen Center for Physics, 
USA, 12-16 January 2004, vol. 328 of ASP Conference Series, 3-18, (Astronomical Society 
of the Pacific, San Francisco, U.S.A., 2005). 

[289] Stairs, I.H., "Prospects of measuring trigonometric parallax for PSR B1534+12", personal 
communication . 

[290] Stairs, I.H., Observations of Binary and Millisecond Pulsars with a Baseband Recording 
System, Ph.D. Thesis, (Princeton University, Princeton, U.S.A., 1998). 

[291] Stairs, I.H., "Testing General Relativity with Pulsar Timing", Living Rev. Relativity, 6, lrr- 
2003-5. (2003). URL (cited on 22 January 2005): 
http : //www . livingreviews . org/lrr-2003-5 

[292] Stairs, I.H., "Pulsars in Binary Systems: Probing Binary Stellar Evolution and General 
Relativity", Science, 304, 547-552, (2004). 

[293] Stairs, I.H., Arzoumanian, Z., Camilo, F., Lyne, A.G., Nice, D.J., Taylor, J.H., Thorsett, 
S.E., and Wolszczan, A., "Measurement of Relativistic Orbital Decay in the PSR B1534+12 
Binary System", Astrophys. J., 505, 352-357, (1998). 

[294] Stairs, I.H., Faulkner, A.J., Lyne, A.G., Kramer, M., Lorimer, D.R., McLaughlin, M.A., 
Manchester, R.N., Hobbs, G.B., Camilo, F., Possenti, A., Burgay, M., D'Amico, N., Freire, 
P.C.C., and Gregory, P.C., "Discovery of three wide-orbit binary pulsars", Astrophys. J., 
submitted, (2005). 

[295] Stairs, I.H., Manchester, R.N., Lyne, A.G., Kaspi, V.M., Camilo, F., Bell, J.F., D'Amico, N., 
Kramer, M., Crawford, F., Morris, D.J., McKay, N.P.F., Lumsden, S.L., Tacconi-Garman, 
L.E., Cannon, R.D., Hambly, N.C., and Wood, P.R., "PSR J1740-3052: a pulsar with a 
massive companion", Mon. Not. R. Astron. Soc, 325, 979-988, (2001). 

[296] Stairs, I.H., Manchester, R.N., Lyne, A.G., Kramer, M., Kaspi, V.M., Camilo, F., and 
D'Amico, N., "The massive binary pulsar J1740-3052", in Bailes, M., Nice, D.J., and 
Thorsett, S.E., eds., Radio Pulsars, Proceedings of meeting held at Mediterranean Ago- 
nomic Institute of Chania, Crete, Greece, 26-29 August 2002, vol. 302 of ASP Conference 
Series, 85-88, (Astronomical Society of the Pacific, San Fransisco, U.S.A., 2003). 

[297] Stairs, I.H., Splaver, E.M., Thorsett, S.E., Nice, D.J., and Taylor, J.H., "A Baseband 
Recorder for Radio Pulsar Observations", Mon. Not. R. Astron. Soc, 314, 459-467, (2000). 
Related online version (cited on 19 March 2001): 
http : //arXiv . org/abs/astro-ph/9912272 

[298] Stairs, I.H., Thorsett, S.E., Taylor, J.H., and Wolszczan, A., "Studies of the Relativistic 
Binary Pulsar PSR B 1534+12: I. Timing Analysis", Astrophys. J., 581, 501-508, (2002). 

[299] Standish, E.M., and Newhall, X.X., "New accuracy levels for solar system ephemerides. (Lec- 
ture)", in Ferraz-Mello, S., Morando, B., and Arlot, J.-E., eds., Dynamics, Ephemerides, and 
Astrometry of the Solar System, Proceedings of the 172nd Symposium of the International 



73 



Astronomical Union, held in Paris, France, 3-8 July, 1995, 29, (Kluwer, Dordrecht, Nether- 
lands; Boston, U.S.A., 1996). 

[300] Stappers, B.W., Bailes, M., Lyne, A.G., Manchester, R.N., D'Amico, N., Tauris, T.M., 
Lorimer, D.R., Johnston, S., and Sandhu, J.S., "Probing the Eclipse Region of a Binary 
Millisecond Pulsar", Astrophys. J. Lett, 465, L119-L122, (1996). 

[301] Stella, L., Priedhorsky, W., and White, N.E., "The discovery of a 685 second orbital period 
from the X-ray source 4U 1820-30 in the globular cluster NGC 6624", Astrophys. J. Lett., 
312, L17-L21, (1987). 

[302] Sternberg Astronomical Institute, Moscow State University, "Scenario Machine Engine", 
project homepage, (2001). URL (cited on 19 March 2001): 
http : //xray . sai .msu.ru/sciwork/ scenario .html 

[303] Stinebring, D.R., Ryba, M.F., Taylor, J.H., and Romani, R.W., "Cosmic Gravitational- Wave 
Background: Limits from Millisecond Pulsar Timing", Phys. Rev. Lett, 65, 285-288, (1990). 

[304] Sudou, H., Iguchi, S., Murata, Y., and Yoshiaki, T., "Orbital motion in the radio galaxy 3C 
66B: evidence for a supermassive black hole binary", Science, 300, 1263-1265, (2003). 

[305] Swinburne University of Technology, "Centre for Astrophysics and Supercomputing" , insti- 
tutional homepage, (2000). URL (cited on 8 December 2000): 
http : // astronomy . swin . edu . au 

[306] Tauris, T.M., and Manchester, R.N., "On the evolution of pulsar beams", Mon. Not. R. 
Astron. Soc., 298, 625-636, (1998). 

[307] Tauris, T.M., and Savonije, G.J., "Formation of millisecond pulsars. I. Evolution of low-mass 
X-ray binaries with P orb > 2 days", Astron. Astrophys., 350, 928-944, (1999). 

[308] Taylor, J.H., "Millisecond Pulsars: Nature's Most Stable Clocks", Proc. IEEE, 79, 1054- 
1062, (1991). 

[309] Taylor, J.H., "Binary Pulsars and Relativistic Gravity", Rev. Mod. Phys., 66, 711-719, 
(1994). 

[310] Taylor, J.H., and Cordes, J.M., "Pulsar Distances and the Galactic Distribution of Free 
Electrons", Astrophys. J., 411, 674-684, (1993). 

[311] Taylor, J.H., and Manchester, R.N., "Galactic distribution and evolution of pulsars", Astro- 
phys. J., 215, 885-896, (1977). 

[312] Taylor, J.H., and Weisberg, J.M., "A New Test of General Relativity: Gravitational Radia- 
tion and the Binary Pulsar PSR 1913+16", Astrophys. J., 253, 908-920, (1982). 

[313] Taylor, J.H., and Weisberg, J.M., "Further experimental tests of relativistic gravity using 
the binary pulsar PSR 1913+16", Astrophys. J., 345, 434-450, (1989). 

[314] The Nobel Foundation, "The Nobel Prize in Physics 1993", institutional homepage, (June, 

2000). URL (cited on 19 March 2001): 

http : //www.nobel . se/physics/laureates/1993 

[315] Thorsett, S.E., Arzoumanian, Z., Camilo, F., and Lyne, A.G., "The Triple Pulsar System 
PSR B1620-26 in M4", Astrophys. J., 523, 763-770, (1999). 



74 



[316] Thorsett, S.E., Arzoumanian, Z., McKinnon, M.M., and Taylor, J.H., "The Masses of Two 
Binary Neutron Star Systems", Astrophys. J. Lett, 405, L29-L32, (1993). 

[317] Thorsett, S.E., Arzoumanian, Z., and Taylor, J.H., "PSR B1620-26: A binary radio pulsar 
with a planetary companion?", Astrophys. J. Lett., 412, L33-L36, (1993). 

[318] Thorsett, S.E., and Chakrabarty, D., "Neutron Star Mass Measurements. I. Radio Pulsars", 
Astrophys. J., 512, 288-299, (1999). 

[319] Thorsett, S.E., and Dewey, R.J., "Pulsar Timing Limits on Very Low Frequency Stochastic 
Gravitational Radiation", Phys. Rev. D, 53, 3468, (1996). 

[320] Toscano, M., Sandhu, J.S., Bailes, M., Manchester, R.N., Britton, M.C., Kulkarni, S.R., 
Anderson, S.B., and Stappers, B.W., "Millisecond pulsar velocities", Mon. Not. R. Astron. 
Soc, 307, 925-933, (1999). 

[321] Trimble, V., "Motions and Structure of the Filamentary Envelope of the Crab Nebula", 
Astrophys. J., 73, 535-547, (1968). 

[322] Tutukov, A.V., and Yungelson, L.R., "The merger rate of neutron star and black hole bina- 
ries", Mon. Not. R. Astron. Soc., 260, 675-678, (1993). 

[323] University of British Columbia, "Ingrid Stairs / UBC Pulsar Group" , personal homepage, 
(2005). URL (cited on 22 January 2005): 
http : //www . astro . ubc . ca/people/stairs 

[324] University of Cagliari, "The Italian Pulsar Group", project homepage, (2005). URL (cited 
on 18 January 2005): 
http : //pulsar . ca. astro . it/~pulsar 

[325] University of California Berkeley Astronomy Department, "Berkeley Pulsar Group" , project 
homepage, (2000). URL (cited on 20 March 2001): 
http : / / astro . berkeley . edu/~mpulsar/ 

[326] University of Manchester, "Jodrell Bank Observatory Pulsar Group", project homepage, 
(1997). URL (cited on 20 March 2001): 
http : //www . jb .man. ac .uk/~pulsar/ 

[327] University of Manchester, "The European Pulsar Network Data Archive" , web interface to 
database, (2005). URL (cited on 17 January 2005): 
http : //www . jb .man. ac .uk/~pulsar/Resources/epn 

[328] van den Bergh, S., and Tammann, G.A., "Galactic and Extragalactic Supernova Rates", 
Annu. Rev. Astron. Astrophys., 29, 363-407, (1991). 

[329] van den Heuvel, E.P.J. , "The binary pulsar PSR J2145-0750: a system originating from a 
ow or intermediate mass X-ray binary with a donor star on the asymptotic giant branch?" , 
Astron. Astrophys., 291, L39-L42, (1994). 

[330] van Kerkwijk, M.H., Bassa, C.G., Jacoby, B.A., and Jonker, P.G., "Optical studies of com- 
panions to millisecond pulsars", in Rasio, F.A., and Stairs, I.H., eds., Binary Radio Pulsars, 
Proceedings of meeting held at the Aspen Center for Physics, USA, 12 Janaury - 16 January 
2004, vol. 328 of ASP Conference Series, 357-369, (Astronomical Society of the Pacific, San 
Francisco, U.S.A., 2005). 



75 



[331] van Kerkwijk, M.H., and Kulkarni, S.R., "Spectroscopy of the white-dwarf companions of 
PSR B 0655+64 and 0820+02", Astrophys. J. Lett., 454, L141-L144, (1995). 



[332] van Kerkwijk, M.H., and Kulkarni, S.R., "A Massive White Dwarf Companion to the Eccen- 
tric Binary Pulsar System PSR B2303+46", Astrophys. J. Lett., 516, L25-L28, (1999). 

[333] van Straten, W., High- Precision Timing and Polarimetry of PSR 30^31-^115, Ph.D. Thesis, 
(Swinburne University of Technology, Swinburne, Australia, 2003). 

[334] Vivekanand, M., and Narayan, R., "A New Look at Pulsar Statistics - Birthrate and Evidence 
for Injection", 3. Astrophys. Astron., 2, 315-337, (1981). 

[335] Weisberg, J.M., "The Galactic Electron Density Distribution" , in Johnston, S., Walker, M.A., 
and Bailcs, M., eds., Pulsars: Problems and Progress (IAU Colloquium 160), Proceedings of 
the 160th Colloquium of the IAU, held at the Research Centre for Theoretical Astrophysics, 
University of Sydney, Australia, 8-12 January 1996, vol. 105 of ASP Conference Series, 
447-454, (Astronomical Society of the Pacific, San Francisco, U.S.A., 1996). 

[336] Weisberg, J.M., and Taylor, J.H., "The Relativistic Binary Pulsar B1913+16: Thirty Years 
of Observations and Analysis", in Rasio, F.A., and Stairs, I.H., eds., Binary Radio Pulsars, 
Proceedings of meeting held at the Aspen Center for Physics, USA, 12 Janaury - 16 January 
2004, vol. 328 of ASP Conference Series, 25-31, (Astronomical Society of the Pacific, San 
Francisco, U.S.A., 2005). 

[337] Wex, N., "The second post-Newtonian motion of compact binary-star systems with spin", 
Class. Quantum Grav., 12, 983-1005, (1995). 

[338] White, N.E., and Zhang, W., "Millisecond X-Ray Pulsars in Low-mass X-ray Binaries", 
Astrophys. J. Lett., 490, L87-L90, (1997). 

[339] Wijnands, R., "Accretion-Driven Millisecond X-ray Pulsars", in Lowry, J. A., ed., Trends in 
Pulsar Research, (Nova Science, Hauppauge, U.S.A., 2005). Related online version (cited on 

23 January 2005): 

http://arXiv.org/abs/astro-ph/0501264. in press. 

[340] Wijnands, R., and van der Klis, M., "A millisecond pulsar in an X-ray binary system", 
Nature, 394, 344-346, (1998). 

[341] Will, CM., "Einstein's Relativity Put to Nature's Test: A Centennial Perspective", lecture 
notes, American Physical Society, (1999). URL (cited on 8 December 2000): 
http : //www . apscenttalks . org/pres_masterpage . cf m?nameID=99 

[342] Will, CM., "The Confrontation between General Relativity and Experiment", Living Rev. 
Relativity. 4. lrr-2001-4, (2001). URL (cited on 24 January 2005): 
http : //www. livingreviews . org/lrr-2003-5 

[343] Wolszczan, A., "A nearby 37.9 ms radio pulsar in a relativistic binary system", Nature, 350, 
688-690, (1991). 

[344] Wolszczan, A., "Confirmation of Earth-mass planets orbiting the millisecond pulsar PSR 
1257+12", Science, 264, 538-542, (1994). 

[345] Wolszczan, A., Doroshenko, O.V., Konacki, M., Kramer, M., Jessner, A., Wielebinski, R., 
Camilo, F., Nice, D.J., and Taylor, J.H., "Timing Observations of Four Millisecond Pulsars 
with the Arecibo and Effelsberg Radio Telescopes", Astrophys. 3., 528, 907-912, (2000). 



76 



[346] Wolszczan, A., and Frail, D.A., "A planetary system around the millisecond pulsar PSR 
1257+12", Nature, 355, 145-147, (1992). 

[347] York, T., Jackson, N., Browne, I.W.A., Wucknitz, O., and Skelton, J.E., "The Hubble con- 
stant from the gravitational lens CLASS B0218+357 using the Advanced Camera for Sur- 
veys", Mon. Not. R. Astron. Soc, 357, 124-134, (2005). 



77 



